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Abstract 



In this thesis we investigate the evolution of open quantum systems in the presence 
of initial correlations with an environment. In the presence of such initial correlations 
the standard formalism of describing evolution by completely positive trace preserving 
(CPTP) quantum operations can fail and non-completely positive (non-CP) maps may 
be observed. 

We investigate a new method of classification of correlations between a system and en- 
vironment using the so-called quantum discord. We found an issue with this classification 
as quantum discord is not a symmetric quantity between exchange of systems. This leads 
to ambiguity in classifications — states which are both quantum and classically correlated 
depending on the order of the two systems. 

Quantum process tomography is investigated with regard to non-CP maps. We ex- 
amine two methods of performing tomography, standard quantum process tomography 
(SQPT) and ancilla assisted process tomography (AAPT). We pay particular attention 
to the effect of state preparation in the presence of initial correlations between the sys- 
tem and environment. We find that in the case of SQPT the preparation procedure can 
influence the complete-positivity of the reconstructed quantum operation. We examine 
a recently proposed method of state preparation by using projective measurements, and 
propose our own protocol that uses a single measurement that is followed by unitary ro- 
tations. In the case of the former the evolution can be non-CP while the later will always 
give rise to a CP map. State preparation in AAPT was found always to give rise to CP 
evolution. In addition we briefly investigate a proposed method of performing bilinear 
process tomography. 

We investigate the effect of statistical noise in process tomography tomography, and 
how it can result in the identification of a map as non-CP when the evolution should 
be CP. We study a probability distribution for quantum operations reconstructed by 
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process tomography in the presence of statistical noise. We found that the variance of 
the distribution for reconstructed processes is inversely proportional to the number of 
copies of a state used to perform tomography. As a result, by increasing this number 
one can distinguish between the distributions of CP and true non-CP processes with a 
high degree of accuracy. Finally, we detail an experiment using currently available linear 
optics quantum computation devices to demonstrate non-CP maps arising in SQPT. 
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Introduction 



One hundred years after its discovery and preliminary development, quantum physics 
continues to offer both new insight into the natural world, and exciting technological 
developments. Techniques exploiting quantum mechanical effects to store, manipulate, 
transmit and process information have been developed in unison with extraordinary sci- 
entific breakthroughs in our understanding of quantum physics. This knowledge and its 
technological applications are encapsulated in Quantum information theory. Quantum 
information is a new branch of research which seeks to develop information processing 
technologies based on quantum mechanical effects. Not only has it found application in 
information centric areas such as totally-secure communication and computation, but it 
has also given remarkable insights into other branches of physics such as relativistic field 
theory and condensed matter physics pQ. 

In this project we shall investigate the theoretical framework which underpins our 
understanding of the processes which occur in quantum information. We will begin by 
briefly introducing the concept of an open quantum system, and how its evolution is 
described. Following this we shall discuss in some detail the fundamental mathematics and 
physics used in quantum information theory to describe the dynamics of quantum systems. 
In particular, we deal with the concept of describing the state of a quantum system by its 
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density matrix, and its evolution by quantum operations using the formalism of completely 
positive maps. We shall also describe the important technique of quantum tomography 
which is used to experimentally measure and determine the state and dynamics of a 
quantum system. Finally we shall discuss the limitations of completely positive maps in 
describing quantum operations and recent developments which attempt to address these 
problems. 

Quantum systems fall into two categories, they are either considered to be open, or 
closed. Closed quantum systems are completely isolated from their environment and hence 
have no interactions with their surroundings. The dynamics of these systems are directly 
described by the postulates of quantum mechanics [2]. More interesting however, are open 
quantum systems. These are systems which exhibit some degree of interaction with their 
environment — all real world systems are of this type to some degree. Understanding 
how these systems evolve is of considerable interest as it is imperative for the construction 
and operation of actual quantum devices to implement quantum information processing. 

In many cases we are not interested in the detailed time-evolution of a process, just 
in the end result. In this case quantum dynamics are described by quantum operations, 
which are maps from a valid initial quantum state to a valid final quantum state for our 
system of interest. According to the standard formalism, quantum operations satisfy two 
requirements: that they are Completely Positive (CP), and Trace Preserving (TP) [3J. 
This is discussed in detail later. Despite being indispensable for much of the theory of 
quantum information science, this formalism has significant limitations. The approach 
based on CP maps assumes that there are no initial correlations between the system and 
its environment. If initial correlations are present the output state of the system after a 
valid physical process may be predicted to have a negative probability [4J. In this case the 
quantum operation is a Non- Completely Positive (non-CP) map. Since a valid physical 
process requires non-negative probabilities, the CP map formalism has failed to model 
the physics involved. 

The problem is that correlations with an environment arise naturally and may be 
present in any experiment. Recent research has focused on trying to classify different 
classes of correlations, however the role of these correlations and a general formalism for 
non-CP maps are not well understood. Unlike the situation with CP maps, only partial 
results are known (5j E] . 



1.1 Overview of the Thesis 
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1.1 Overview of the Thesis 

In this thesis we investigate certain situations where non-CP maps arise from initial 
correlations between an open system and the environment. The layout of material covered 
is as follows. 

In Chapters [2] and [3] we introduce the essential mathematical background required 
for the remaining chapters of this thesis. We introduce the concept of quantum infor- 
mation and encoding information into the state of a quantum system. We describe the 
representation of the state of a quantum systems by its density matrix, and some useful 
decompositions of this representation. We introduce the standard formalism for describing 
the evolution of open quantum systems by quantum operations, represented by completely 
positive trace preserving (CPTP) maps. In particular we describe several different math- 
ematical representations for CPTP that were encountered in the literature. We provide 
clear summary of the relationships between these representations as this was not found 
in the literature. 

In Chapter [4] we introduce non-CP maps and discuss how they can be described using 
assignment maps [5J. We also investigate the classification of initial correlations between a 
system and environment. We follow a recently proposed method of classifying correlations 
based on a quantity called quantum discord [SI El IE] ■ However, we found problems with 
this approach as quantum discord is not a symmetric quantity. We show this by counter- 
example. The asymmetry prompted us to conjecture that a theorem [6] concerning the 
relationship between CP evolution and quantum discord does not hold when we exchange 
the roles of two systems. 

In Chapter [5] we investigate the quantum process tomography, which is the process 
of characterizing an unknown quantum operation. We introduce two methods, standard 
quantum process tomography (SQPT) and ancilla assisted process tomography (AAPT). 
As a preliminary step we describe the characterization of unknown quantum states by 
state tomography. Here we interpret these schemes in a notation consistent with that 
used in Chapter [3} We also discuss sources of statistical noise in process tomography, 
how this can give non-CP results, and current schemes for dealing with the the noise by 
maximum likelihood process tomography. We raise the issue of this scheme not being able 
to distinguish between true non-CP results arising from initial correlations, and non-CP 
results due to statistical noise. 

In Chapter[6]we investigate the effect of initial correlations on several state preparation 
techniques used in process tomography. We begin by investigating preparation schemes 
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proposed in the literature, and then we propose a new scheme. We extend the previous 
work by relating state preparation to non-CP evolution by contrasting the effect of the 
different schemes when we take into account initial correlations. At the end of the chapter 
we introduce a bilinear process tomography scheme proposed in Ref [H] . However, we find 
one of the asserted properties of the matrix describing bilinear evolution incorrect. We 
provide the corrected result. 

In Chapters [7] we investigate the statistical noise, introduced in Section 5.3 arising 



in optical implementations of process tomography. We propose original ideas for distin- 
guishing between non-CP results arising from from noise, and those that legitimately arise 
from initial correlations. 

In Chapter [8] we bring together our results from the preceding chapters to outline 
an original experiment which could be performed with currently available linear optical 
techniques to demonstrate how the different state preparation procedures, introduced in 
Sections 6.1.2 and 6.1.3[ can result in non-CP evolution. 

Finally, in Chapter [9] we review the main results from this thesis and discuss directions 
for future investigations which have arisen from our work. 



2 

Quantum Information 



We begin with some essential background in the field of quantum information science. 
Before we can describe the evolution of open quantum systems in the presence of initial 
correlations, we must first introduce some fundamental tools used in quantum information 
theory. In this chapter we will briefly introduce the idea of storing information in the 
state of a quantum system, and the basic mathematics required for the description of 
these quantum systems and their states. 

2.1 Qubits 

In classical information theory all information is described by strings of bits, each of 
which can have a value of or 1. Computation is then achieved by performing logical 
operations on these strings. In quantum information theory one represents information as 
quantum bits (qubits). A qubit represents a two-level quantum system whose state space 
is spanned by the state vectors |0) and |1)R Actual two-level physical systems which make 

1 Strictly speaking, a qubit is represented mathematically by a two-dimensional complex Hilbert Space 
with orthonormal basis {|0) , |1)}. This choice of basis vectors is known as the computational basis , and 
the vectors are labeled to correspond to their classical bit counterparts. 
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good candidates for a qubit include the spin of a particle (up or down), the polarization 
of a photon (vertical or horizontal) , or the state of a two level atom (ground or excited). 

An important distinction between a qubit and a classical bit is that unlike its classical 
counterpart, the state of a qubit can be any normalized complex linear combination of its 
basis states: 

|^)=a|0)+/9|l), \a\ 2 + \(3\ 2 = 1. (2.1) 

Physically, if we were to measure the qubit with a device capable of detecting either 
|0) or |1), then the result would be |0) with probability |a| 2 , or |1) with probability |/3| 2 . 
This is known as Bom's Rule and it connects the complex coefficients a and (3, known as 
amplitudes, with the probability of measuring results. 

To model larger systems consisting of multiple qubits we introduce further notation. If 
we have two qubits labeled A and B with basis states |0) A , and |0) B , |1) B respectively, 
the basis states for the composite system are given by the tensor product of those for the 
individual systems. Hence, for example, in a 2-qubit system the basis states are 



{|00> , |01> , 1 10) , |11)}, where |00) = |0>, ® |0) 



B ■ 



this will be discussed in more detail in Section 2.6 We will now introduce several impor 



tant mathematical tools needed for the description of the states of a quantum system. 



2.2 The Density Matrix 

To this point we have described our quantum systems in the language of state vectors. 
However, there is a more general approach using what is known as the density operator or 
density matri^ The density matrix approach is mathematically richer than the state vec- 
tor approach, and provides a more convenient framework for many scenarios encountered 
in quantum information science. 

To be precise, if a quantum system is in one of several states \ipi) with corresponding 
probability pi, the density matrix, p, for the system is defined by: 

i 

2 When dealing with finite-dimensional quantum systems the terms operator and matrix can be used 
interchangeably. This is because given a basis for the state space of the system, all linear operators acting 
on the space have a matrix representation with respect to the given basis. 



2.2 The Density Matrix 



7 



If a quantum system is known to be in a state \ip) with unit probability, its density 
matrix is given by p = and it is said to be in a pure state. This case is equivalent 



to the state vector description in Eqn. (2.1). A density matrix defined as a mixture of 



different pure states, as in Eqn. (2.2), is said to be in a mixed state. Note that if the 
density matrix for a qubit is given by p = 1/2, where I is the identity matrix, p is said to 
be in a maximally mixed state. 

For example, if we had a single qubit in the state 

P = Po |0)(0| +pi |1)(1| , where p +p 1 = l, 

then the probability of measuring a |0) or |1) is po or p\ respectively. If po = Pi = 1/2, 
this is a maximally mixed state, while if either po = or p\ = this will be a pure state. 
A measure of the purity of a state p is given by Tr(p 2 ). Pure states satisfy Tr(p 2 ) = 1 
while for mixed states Tr(p 2 ) < 1 

For p to be a valid density matrix for our system it is necessary and sufficient that it 
satisfies two conditions. 

1. p is a positive-semidefinite matrix: A matrix p is positive-semidefinite if it is 
positive and hermitian. 

• A matrix p is hermitian if pt = p, where t denotes conjugate-transposition, 
(p^)mn = Pnm- This is equivalent to requiring that all the eigenvalues of p are 
real. 

• A positive matrix p satisfies \ p \ip) > for any vector This is equivalent 
to requiring that all the eigenvalues of p are non-negative. 

2. p has unit trace: That is to say Tr(p) = p mm = 1. Along with the requirement 
of positive-semidefinite, this is equivalent to the eigenvalues of p being real, non- 
negative and summing to 1. 

The requirement of positivity is related to measurement probabilities. If is a pure 
state, then (tp\ p \ip) is interpreted as the probability of measuring |?/>) given p. This is why 
we require it to be non-negative. Requiring Tr(p) = 1 guarantees that the probabilities 
add up to 1. In general we will refer to a positive-semidefinite matrix simply as a positive 
matrix, denoted by p > 0. We also note that positivity implies hermiticity in a complex 
Hilbert space. 
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2.3 Bloch Sphere 

We can visually represent the state of a qubit as a point on the Bloch sphere. An example 



is shown in Fig. (2.1). For convenience, the axes of the sphere are labeled to correspond 



to polarization states of a photon. The correspondence between the polarization states 



and the computation basis {|0) , |1)} is shown in Table (2.1). 




|D> 



FIGURE 2.1: The pure state \D) = (|0) + |1))/V2 represented on the Bloch sphere 



Photon polarization 




State vector in 


state vector 




computation basis 


Horizeontal 


\H) = 


|o> 


Vertical 


\V) = 


|1) 


Diagonal 


\D) = 


I+> = ^(|0> + |1» 


Anti-Diagonal 


\A) = 


l-> = ^(|o>-|i» 


Right Circular 


\R) = 


I+<> = ^(|0>+<|1» 


Left Circular 


\L) = 


H> = ^(|o>-<|i» 



Table 2.1: Correspondence between polarization state of a photon and the computational 
basis for a two-level quantum system 



Any pure state vector of a qubit can be parameterized by two angles 9, 4> by 

|V) = cos0|O) +e i0 sin0|l) . 



2.4 Spectral Decomposition 
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This corresponds to a point on the surface of the Bloch sphere specified by the vector 
a e IR 3 , where a = (01,02,^3) = (sin 26 cos </>, sin 26 sin <p, cos 26). The vector a is called 
the Bloch vector for the state \(f>). 

A density matrix can also be represented in terms of its Bloch vector, a, except now 
a need not be of unit length. Pure states are represented by points on the surface of the 
sphere and have \a\ = 1. Mixed states are represented by points inside the Bloch sphere 
and have \a\ < 1. The maximally mixed state corresponds to the center of the Bloch 
sphere and has a — 0. For an arbitrary density matrix p, the Bloch representation is 
given by 

1 3 



8=1 



where <r, are the Pauli matrices 



o~i = X 



1 

1 



a 2 = Y 



-i 

1 



0-3 



Z 



1 
-1 



(2.3) 



The set {I, <7i, 02, 03} is an orthogonal basis for all 2 x 2 complex matrices with respect 
to the Hilbert- Schmidt inner product (A, B) = Ti(A^B). We can recover the elements of 
the Bloch vector for a state p by «j = Tr(po"j). 



2.4 Spectral Decomposition 

A particularly useful result from linear algebra is the spectral theorem [TU] . This theorem 
applies to normal matrices. A matrix M is normal if and only if M^M = MM^ . Hermitian 
and unitary matrices are both subsets of normal matrices. 

The spectral theorem states that for any normal matrix M on a vector space V there 
exists an orthonormal basis for V which diagonalizes M. The basis is given by the nor- 
malized eigenvectors, |e;), of M, and the entries of the resulting diagonal matrix Dm are 
the eigenvalues, Aj, of M. Hence we can form the spectral decomposition of M, 

M = E M D m E M = z ~2\ \ei)(ei\ , 

i 

where | e^) are the columns of E M . 

In the case of density matrices, the eigenvalues are real, non-negative and sum to 1. 
This is an important result as it means that every density matrix can be expressed as a 
convex-linear sum of its eigenstates | e^) (e^ | . A sum p = J2jPjPj * s convex-linear if the 
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coefficients pj are all non- negative and sum to 1. Conversely, a density matrix is a pure 
state if it only has one non-zero eigenvalue. 

2.5 Postulates Of Quantum Mechanics 

Quantum mechanics can be formulated using density matrices and the following four 
postulates PQ: 

• Postulate 1: The state of any closed physical system is completely described by 
its density matrix. The density matrix, p, is a positive matrix which acts on the 
state space of the system and has trace one. Here a closed system refers to a system 
which does not interact with its environment. If the system is in a state pi with 
probability pi, the systems density matrix is given by ^ZiPiPi- 

• Postulate 2: Closed quantum systems evolve under unitary transformation. That 
is, if a system evolves from p to p' at times t\ and ti respectively, the states are 
related by a unitary operator U = U(t2,ti) by the equation 

p? = UpU\ (2.4) 

where U* is the adjoint of the operator U (an operator U is unitary if and only 
if WU = I). This is equivalent to saying that closed quantum systems evolve 
according to Schrodinger's equation. 

• Postulate 3: Quantum measurements are described by a collection of positive op- 
erators, {M m }, acting on the state space of the system being measured. These 
operators satisfy the completeness relation M^Mm = I. If, prior to measure- 
ment the system is in state p, the probability that result m occurs is given by 

p(m) = Tr{MlM m p), 

and the post measurement state of the system is 

M mP Mj 
Pm Ti{MlM m p) ' 

• Postulate 4: For a composite physical system, the state space is given by the tensor 
product of state spaces of its component systems. If we have n constituent systems, 
with the i system being prepared in state pi, the joint state of the composite system 
is pi ® p 2 (8) . . . <g> p n . 



2.6 The Tensor Product 
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The tensor product mentioned in Postulate 4 is a mathematical operation for combining 
two or more matrices of arbitrary dimension into a larger block matrix [TT]. If we consider 
two matrices A and B, with matrix elements a mn , and bki respectively, then their tensor 
product is given by the block matrix 

^ a\\B ai 2 B . . . a\ n B \ 

ao\B ' ■ ■ : 

A®B = 

\ a mlB ■ ■ ■ dmnB ) 

This operation can be applied to both density matrices and state vectors, we treat the 
latter as d X 1 matrices. 

For matrices (or vectors) A, B, C, D and scalar a the tensor product satisfies the 
following relations: 



Left and right distributivity: 
Bilinearity: 
Associativity: 
Multiplication: 
Trace distributivity: 
Adjoint, Transpose, Inverse, Complex-conjugate: 



(A + B) ® (C + D) = A® B + C ® D 
(aA) ® B = A® (aB) = a(A <g> B) 
{A®B)®C = A®{B®C) 
{A ®B)(C®D) = {AC) ® (BD) 
Tr(A ® B) = Tr(A)Tv(B) 
(A ® B)° = A°(g)B° where ° 



_t T -1 

5 ) 



2.7 Composite Quantum Systems 

We shall now introduce notation for describing composite systems of two qubits. However, 
by repeated application one can extend this to systems of many qubits. The Bloch 
representation of a single qubit can be generalized to a system of two qubits, A and B, 
by taking the tensor products of the basis elements {I, o~±, 02, (73}. In this notation we can 
express an arbitrary two qubit state as 

PAB = 4 Yl ^ ® 1 + ai0i ® 1 + Pi 1 ® °i + 7v ai a ^ ■ ( 2 - 5 ) 

ij 

where a = (oti, a 2 , 03) and (3 = fi 2 , P3) are the Bloch vectors for the reduced states of 
system A and B respectively, and 7^ are real parameters describing correlations between 
the systems. 
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If the state of the combined system can be written as pab = Pa® Pb, then pab is said 
to be simply separable, or a product state. If pab can be expressed as a convex-linear sum 
of product states, pab = YliPiPt ® pf ■> then p AB is said to be separable. 

An important feature distinguishing quantum systems from classical ones is that when 
we combine two systems A and B, arbitrary states of the joint system AB may have 
correlations present which cannot be described classically. In the cases where pab is not 
separable, the state is said to be entangled. Entanglement is a feature unique to quantum 
physics and is the reason for many of the intriguing results of quantum information, such 
as quantum teleportation and totally secure communication [1] . We consider correlations 
in more detail in Chapter |4| 



2.8 Partial Trace 

If we know the state of a composite system, the state of one of the component systems 
can be described using the partial trace. Consider two systems, labeled A and B, with 
the composite system described by the density matrix pab- The density matrix for the 
system A is then given by the reduced density matrix 

p A = Ttb(pab)- (2.6) 

Here Tr# is called the partial trace over system B. 

If we have a system of two arbitrary matrices A ® B, then Trg is defined as: 



Tyb(A® B) = ATr(B). 



(2.7) 



By linearity the above expression can be extended to a more general correlated states 



Pab, which we can express in the form of Eqn. (2.5 ) or its higher dimension generalization. 



From here we can directly apply the partial trace as defined by Eqn. (2.7). Since the Pauli 
matrices are traceless the reduced states are 



Pa 



OLiOi 



PB 



jO~i I , 



1=1 



which is the Bloch sphere representation from Section 2.3 

A second method of calculating the partial trace of a density matrix pab is in terms 



2.9 Schmidt Decomposition 
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of any orthonormal basis {\i)} for the system we wish to trace over. In this case we have 



p A = TtbPab 

= ^2 (if Pab 

i 

= ^2 Tl (pAB\i)(i 

i 

AB )mi,ni ■> 



(2.8) 
(2.9) 

(2.10) 

(2.11) 



where {\i)} is an orthonormal basis for B and the superscript B denotes that the operator 
acts only on system B, i.e. if B is a qubit 



\i) B = 1® \i) 



Ui 





Vi 








Ui 





Vi 



where \i) 



Ui 
Vt 



Similarly we can define operators \i) = \i) ® I to trace out system A. 



2.9 Schmidt Decomposition 

Another important result which we will briefly mention is the Schmidt decomposition pp. 
This theorem states that for any pure state of a composite system AB, there exist 
orthonormal bases {|^a)} and for systems A and B respectively such that 

H) = y^A^ \i A ) ® \i B ) . 

i 

The coefficients Aj, called the Schmidt coefficients, are real, non-negative numbers satisfy- 
ing Y^i A? = 1. The number of non-zero Schmidt coefficients is called the Schmidt number 
for and represented by Sch{\ijj)). 

The Schmidt number can be used as a measure of entanglement. For any state \ifj), 
Sch(\ifj)) = 1 if and only if is a product state, while Sch(\%p)) > 1 if and only if \ip) is 
entangled. 

Generalizing this to density matrices gives 

P = Va){3a\ ® |«b>0'bI • 

i,3 
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The density matrices for the reduced states are given by pa = Ttb(p) = X^Af 

and pB = Tta(p) = $^Af Hence the reduced states pa and ps have the same 

eigenvalues. For a density matrix a Schmidt number of 1 corresponds to a simply separable 

state. 



2.10 Vectorization of Matrices 

Since the space of matrices is a linear space, we can represent density matrices as vectors 
on a higher dimensional Hilbert space know as a Liouville space [11]. This is useful for 
describing certain topics encountered later on. This process is called vectorization, and 
transforms an n x m matrix M into a (mn) x 1 column vector denoted by vec(M) or 
\M). This is done by stacking the columns of M to form a vector, so if M^ = m^ then 
vec(M) = |M) = [m 11 ,...,m ml ,m 1 2,...,m m2 ...m ln ,...m mn } T . 
In the case of a single qubit we have 



P 



Pll Pl2 
P21 P22 



then \p) 



Pn 
Vix 

Pl2 
P22 



The inverse of vectorization is the mat function, so mat \p) = p. We will now list some 
useful properties of vectorized matrices [IT]. Let A, B and C be complex square matrices. 
Then, 

1. The vector inner product is the Hilbert-Schmidt inner product: (A\B) = Tr(A^B). 

2. The Schmidt number of \A) is equal to the rank of A. 

3. Vectorization is linear: For matrices Aj and scalars a iy |^ ajAj) = ^\ a, 

4. The action of a matrix A <g> I on a vectorized matrix |S) is A <g> I |S) = |5y4 T ) 

5. More generally, we have the identity \ABC) = (C T A) \B). 

6. For computational basis if A = a^ then \A) = a^- |j) ® 

7. Finally, for any matrix A, \A) = A T ®I\cp) — I ® A \cp) where \cp) = J2i I*) ® K)- 



3 

Quantum Operations 



Now that we have introduced the essential mathematical tools for describing the states 
of quantum systems, we move onto the description of their evolution. Postulate 2 is 
only sufficient for describing the dynamics of closed quantum systems. However, in the 
real world there will always be some degree of interaction between a system and the 
environment. We require a more robust framework to describe the dynamics of these 
open quantum systems. The standard mathematical formalism for the evolution of open 
systems is quantum operations, though this formalism does have its limitations which we 
will discuss. 

A quantum operation is a map, 8, acting on the state space of a quantum system. 
The map describes how quantum states, represented by density matrices p, transform. 



A formal definition will be given in Section 3A_ Two elementary examples of quan- 
tum operations are unitary transformation A(p) = UpW, and quantum measurement 
A m (p) = MrnpM^/TrlM^ M p], which were mentioned in postulates 2 and 3, respectively, 
of Section 12.51 

There are several different approaches for describing general quantum operations which 
we will introduce in the following sections. The formal definition of a completely-positive 
trace-preserving map (CPTP) is based on a set of physically motivated axioms such 
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an operation should satisfy. There are then several useful mathematical representations 
for describing CPTP maps, these include the Kraus representation, process matrix, and 
the superoperator. All these methods are equivalent under the condition of complete 
positivity, however there are certain advantages to each approach [Tj, [3J [TJl H3] • Finally 
there is a physical interpretation where we consider the system and environment together 
to be a closed system. 



3.1 Completely Positive Maps 

First we shall introduce the formal definition for a quantum operation. This definition 
is based on a set of axioms encompassing the physical constraints the evolution of a 
quantum system should satisfy. It mostly follows from the properties of density matrices 



in Section 2^2, which a map describing quantum evolution should preserve. If the output 
state is to be a valid density matrix we would expect the map to be trace preserving, 
convex-linear, and positive. It turns out that requiring the map to be positive is not strong 
enough, and instead complete positivity is required. The reason for this requirement is 
explained below. 

A quantum operation for a system S is a map, 

S:p^S(p), (3.1) 

acting on density matrices for a system S, which satisfies the following axioms: 

1. £ is trace preserving. That is to say Tr[£(p)] = Tr(p). 

2. £ is a convex-linear map. That is 

^ I ^PiPi I = y^Pi£{Pi) where p { > OVi, and = 1 (3.2) 

V i / i i 

3. £ is a completely-positive map. This means that £ is positive (maps positive matrices 
to positive matrices), and also if we introduce an auxiliary system of arbitrary 
dimension then the map X ® £ on the joint system is positive, where X is the 
identity map on the auxiliary system. 

The trace preserving property of £ ensures that the output density matrix will satisfy 
the trace condition. Convex-linearity ensures that the transformation of a mixed state 
is equivalent to the probabilistic sum of the transformations of constituent pure states. 



3.2 Mathematical Representations of Completely Positive Maps 
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Finally, complete positivity of the map ensures that the positivity condition is satisfied by 
the output state of any combined system. Requiring that £ be only a positive map instead 
of completely positive is not sufficient as it does not guarantee that the transformation of 
a positive state on a composite system is itself positive. We will show this by example. 

Consider the transposition map T(p) = p T . Since for any density matrix p T > 0, T is 
a positive map. However if we consider the action of X ® T on the maximally entangled 
2-qubit state -4= ( 1 00) + |11)), the density matrix for the final state of the joint system is 





/ 1 








o\ 


1 








1 





2 





1 










U 








1 / 



which has a negative eigenvalue of —1/2 and so is not a physical density matrix for the 
joint system. 

A map satisfying these three axioms is referred to as a completely positive, trace 
preserving (CPTP) map. It is possible to relax the trace condition to Tr [£(/?)] < Tr[p], 
allowing for completely positive trace decreasing maps. We will only be concerned with 
CPTP maps and for the remainder of this thesis, except when explicitly specified, we shall 
refer to maps satisfying all three of these axioms simply as completely-positive maps (CP 
maps) . 



3.2 Mathematical Representations of Completely Pos- 
itive Maps 

Now that we have defined quantum operations to be completely positive maps, we need 
a mathematical representation for them. There are several different representations used 
in literature, and we will now introduce three of them. These are the Kraus represen- 
tation [3], the process matrix [12], and the superoperator [13]. We will also reveal the 
relationships between these representations, and how each can be converted into the oth- 
ers. Following this, we will show how these representations can arise from considering 
unitary evolution of the system and environment when taken from the point of view of 
the system alone. 
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3.2.1 Kraus Representation 

Our first representation, and the most widely used, is the Kraus representation, and is 
the result of the following theorem [12] . 

Theorem 3.1 (Kraus Representation) A map £ acting on density matrices of a sys- 
tem S is CPTP if and only if there exists a set of operators {K n } acting on the state 
space of system S such that 

£{ P ) = Y,K n pKi where J2 K t K - = L ( 3 - 3 ) 

n n 

This form of expressing a quantum operation is also refereed to as the operator-sum 
formalism, and proof of this theorem can be found in most textbooks concerning quantum 
information (TJ [121 02] • The operators K n are called Kraus matrices and they satisfy 
Y2 n K\K n = I, which is known as the completeness relation. 

An elementary example is unitary evolution, where we only have one Kraus matrix. In 
general the evolution described here need not be unitary. The Kraus representation allows 
us to completely characterize the dynamics of an open system S by a map acting only on 
the state space of S. In other words we do not have to explicitly consider properties of 
environment, they are all accounted for in the Kraus operators {K n }. Another feature of 
the Kraus representation is that it is not unique. This can be useful as different system- 
environment interactions may still give rise to the same reduced dynamics on the system. 

3.2.2 Process Matrix 

Our second method describing a completely positive map £ acting on a system S is with 
the process matrix, A^. We define the unique process matrix for a CPTP map £ as a 
matrix on a fictitious space A® S, where the auxiliary system (ancilla) A is a copy of S. 
It is given by 

A e = y %2E ij ®e(E ij ), (3.4) 

where E^ = is the matrix with a 1 in the (ij) th entry and zeros elsewhere. 

The evolution of a state p by a quantum operation £ can then expressed in terms of 
the process matrix by 



£{p) = Tr A [A s {p T ® I)] = Tr A [{p T <g> I)A £ ] 



(3.5) 



3.2 Mathematical Representations of Completely Positive Maps 
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The proof of this is given in Appendix A.l.l| p 



There is an equivalence between the completely positivity of a quantum operation £ 
and properties of its process matrix. This equivalence is given by the following theorem. 

Theorem 3.2 (Process Matrix) A quantum operation £ on a d- dimensional system S 



is CPTP if and only if its process matrix As, defined by Eqn. (3.4), satisfies 

1. As is positive- semidefinite (As > 0). 

2. Tr B (A £ )=I. 

For a proof we refer the reader to [12J or [13J. 



From Eqn. (3.4) it is easy to see how to construct a process matrix from a Kraus 
representation. If we wish to switch representations in the other direction we note that 
for a completely positive map £ the process matrix is positive and hermitian. Hence 
As has non-negative eigenvalues and orthonormal eigenvectors, \ and | e^) respectively. 
We construct the Kaus matrices by 



K n = v A n mat | e^) . 
The proof of this can be found in Appendix |A.1.2 



Process Matrices as States 

As has all the properties of a density matrix except the normalization |^} It is a positive 
hermitian matrix, however it is d 2 x d 2 instead of d x d, and has trace d instead of 1. This 
similarity is due to a relationship known as the Jamiolkowski isomorphism |12| . which 
gives a direct correspondence between process matrices for a rf-dimensional system S, and 
density matrices for a <i 2 -dimensional system A <g> S, where we have introduced an ancilla 
A of the same dimension as S. 
To see this we note that 



i,j hi 

where \<p) = ® K))/ 1 /^ is a maximally entangled state of the joint system A <g> S. 



1 One could also define the process matrix by A^ = J^i j ® Eij- O ur evolution would then be 

given by £(p) = Txb[A£(I(S> p T )]- Since we can treat Ag as a bipartite matrix, this change in definition is 
just a swap of the two components A «-> B. The relationship between Af and the other representations 
would also change with respect to this swap. 

2 The normalization is an aesthetic choice in order to simplify the related equations. 
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Figure 3.1: The relationship between a process matrix and a density matrix by the Jami- 
olkowski isomorphism 



Thus the correspondence is given by p £ = (I ® £)(\<f>){</>\) = A/d. This can be 
thought of as introducing an ancilla system A as illustrated in Fig. (3.1). One important 
application of this approach is in process tomography which we investigate in Chapter [5j 



3.2.3 Linear Superoperator 

Our final mathematical representation of a quantum operator £ is as a linear superoperator 
$£, which acts on vectorized density matrices. This evolution is given by 

\£(p)) = $ £ \p) . (3.6) 

We give a novel definition for the superoperator, which is not found in the literature, 

* e = ^2\E{E ij )){E ij \, (3.7) 



where = is the matrix with a 1 in the (ij) entry and zeros elsewhere. Proof of 



this can be found in Appendix A. 1.3 



While a linear superoperator is a convenient representation of a CPTP map £ , proving 
an arbitrary map is CP directly from $£ is not very illuminating. What is usually done [1_3] 
is to transform $£ into the process matrix A^. This is done by a process called reshuffling. 
Reshuffling is a rearrangement of the matrix elements of the superoperator. It is given by 
{^e)mn,iiv — {^s)^,m,vn- Full details of this procedure can be found in [TH] . 

However, we noticed a more intuitively picture of reshuffling by using our definition for 



the superoperator. If we compare the definitions of $g in Eqn. (3.7) and in Eqn (3.4) 
we can see that this process swaps 



3.3 System-Environment Model 
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The superoperator can easily be calculated from the Kraus representation by vector- 



izing and using the identities in Section |2.10 

\£{p)) 



vec^2K n pKl\ (3.8) 
^K*®K n \p) (3.9) 



n 

so $ = Y2n-^n ® ^n- The reverse process is not so simple. One approach is to reshuffle 
$g to form As, and the apply the procedure for constructing a set of Kraus operators 
from a process matrix. 

3.2.4 Summary of Relationships Between Mathematical Repre- 
sentations 

We have now introduced several different mathematical descriptions for CPTP quantum 
processes. In this section we will briefly summarize the relationships between them. To 
our knowledge such a compact summary is not presented anywhere in the literature though 
it is a useful computational aid. 

For a CPTP quantum operation £ acting on system A there exists a process matrix 



Ag, a superoperator $£, and Kraus operators {K n } as defined in Sections 3.2.2, 3.2.2 



and 3.2.3 respectively. Evolution of density matrix p of system A is then given by 



Tr A [A £ (p T ®I 
S(p) = { mat($ £ \p)), 

H n K n pKl 

To switch between representations we can use the relationships summarized in Ta- 



ble (3.1) 



3.3 System- Environment Model 

Our final approach for modeling open quantum systems is built on an intuitive picture of 
the joint system comprised of our system of interest, and the environment it is interacting 
with. We will denote these by quantum systems S and E, respectively. We also assume 
that the state space of environment is of the same dimension as the system pp. The 
composite system SE is assumed to be a closed system and thus its dynamics are described 



by the postulates of quantum mechanics mentioned in Section 2.5 
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To \ From 


Process Matrix (A^) 


Superoperator 


Kraus Representation (-ftT n ) 


Af = 


A £ 


Reshuffle 

( Ag J ■xnn,\w fim.vn 


^« ® (En K nEijKl) 


$ £ = 


Reshuffled) 






K n = 




Reshuffle^) 
y/Xntnat 





Table 3.1: Relationships between mathematical representations of CP maps. 



Suppose our system is in a state p. If we send this state into a black box n where it 
undergoes evolution while interacting with the environment, the final state S(p) will in 
general not be related to the initial state p by a unitary transformation. To describe the 
evolution of p we must consider the unitary transformation of the whole system. We can 
then recover the transformed state of the system by performing the partial trace over the 



environment. This is illustrated in Fig. (3.2). 



P 



T 



S 



E 




Hi' 



Figure 3.2: Model of an open quantum system as a system interacting with its environ- 
ment m 

To do this we first suppose the environment is initially in a state tq. We also assume 
that the there are no correlations between the system and the environment, and hence 
the state of the combined system is given by the p®r . This is a crucial assumption. The 
final state of the system is then given by: 

£(p)=Tr E [U(j>®T )Ut]. (3.10) 

This leads us to the following theorem 

3 A black box is a device which does something to our system, and the specifics are irrelevant. All we 
are concerned with is is what we put in one side, and what comes out the other. 
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Theorem 3.3 For an open quantum system S with environment E, if the initial state 
on the joint system SE is a product state, then the quantum operation S, as described in 



Eqn. (3.10), will always always be CPTP. 



The proof of this is a shown in Appendix A. 1.4 , and is a matter of constructing a Kraus 
representation for the evolution. 

While this method provides a nice physical picture for the evolution of an open system 
it is necessary to know details of the environment and the systems interaction with it. 
Another limitation of this approach is that we have assumed that the composite system is 
in a product state. This generally requires the experimenter to perform a preparation pro- 
cedure to remove any correlations between the system and environment. The implications 
of this are discussed in Chapter [6| 

If we now consider a more general situation where there are initial correlations between 
the system and environment, which is to say the initial state of the joint system may be 
entangled, or the state of the environment somehow depends on the state of the system, 
then the formalism of quantum operations as CP maps breaks down jU |5j EJ El HHJ HHJ 
rm [18| [19] . In these cases one cannot construct a Kraus representation, and the process 
matrix will have negative eigenvalues. There are no general mathematical representations 
of non-completely positive maps [5]. 

This completes our review of the mathematics used in the standard formalism for 
describing the evolution of open quantum systems. In the next chapter we consider the 
nature of initial correlations between a system and environment, which can cause the 
presented formalism to break down. 
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Non-Completely Positive Maps 



In the previous two chapters we introduced several mathematical tools needed to describe 
the state and evolution of open quantum systems. Open quantum systems were assumed 
to interact with an environment, and the system together with its environment was treated 
as a closed system. The most general description of the state of a system was given by 
its density matrix, and evolution was described by a completely positive trace preserving 
map which acts on the density matrices. We mentioned that the CP formalism can break 
down in the presence of initial correlations between an open system and its environment. 
In this chapter we investigate such correlations between two quantum systems, and their 
effect on quantum evolution. 



4.1 Separable and Entangled States 



As was introduced in Section 3.3, we will be considering S which is interacting with an 
environment system E of equal dimension. We will specifically deal with the case of qubits 
(d — 2), though these results can be generalized to higher dimensional systems. Recall 
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from Section 2/7 that any state of the joint system SE can be written in the form 

^33 

PSE = ^ 2J (I ® I + <8> I + <8> CTj + 7y<7i <g> (Tj) , 

4=1 i=l 

where cr, are the Pauli matrices. The states of the reduced systems are given by 



P8 = \{^-+YL aia ^j ' and p e = \(^ + Yl fajj ■ 

We have mentioned three classes for categorizing the correlations of joint state pse- 
These are 

1. Simply separable: pab = Pa® Pb- 

2. Separable: pab = J2iPiPt ® pf ■ Here pab is a convex-linear sum of simply 
separable states. 

3. Entangled: pab is not separable. 

Traditionally, it was thought that classical correlations (where the system and en- 
vironment are separable) give rise to CP dynamics, while quantum correlations (where 
system and environment are entangled) could give rise to non-CP dynamics [20]. How- 
ever, while simply separable systems always have CP dynamics, it has been shown that 
non-entangled systems can still give rise to non-CP dynamics [3 Ej. We will illustrate 
this with an example. 

4.1.1 Example of a Non-Completely Positive Map 

We will use an example from [5], where the state of a joint system SE is dependent only 
on the reduced state of system S. The initial state of a two qubit system is 

Pse = T I 1 se + ^2 ai(Ti ® 1 + a ^ ai ® ai ' ^ 4 - 1 ) 



4=1 i=l 



where a is a fixed parameter. In this case pse > if < \a\ < | 'A — 3 1 ct | 2 — lj, so 
if we fix the value of a, only certain reduced states of S will give a physical state of the 
joint system. In addition the state pse is always separable for a > 0. This is a case of an 
assignment map which we will formally introduce momentarily. 



4.2 Describing Non-CP Dynamics by Assignment Maps 
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If we evolve this state by a two-qubit unitary rotation 
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This state and unitary evolution lead to non-CP evolution on the first system [5]. This 
can be seen by the eigenvalues of the process matrix taking negative values for certain 



values of 9. The eigenvalues for a fixed value of a = 0.2 is shown in Fig. (4.1). 

A 
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-0.02 



-0.03 







Figure 4.1: Three of the eigenvalues of the process matrix for evolution of pab by U{9). 
Reproduced from [5] with permission. 



This indicates that using entanglement as the definition of quantum correlations is 
too restrictive. A more recent formalism for describing quantum correlations in terms of 



what is known as quantum discord will be discussed in Section 4^3 However, first we 
must describe assignment maps and how they can give a physical meaning to non-CP 
dynamics. 



4.2 Describing Non-CP Dynamics by Assignment Maps 

One approach to characterize quantum operations, CP and non-CP, is with the notion 
of an assignment [2T] or extension [22] map. Let us consider an open quantum system 
S interacting with an environment system E. In the approach used in Section |3.3| we 
started with a state of the joint system p SE with the property ps = Tte(pse), and 
considered its unitary evolution. Any correlations between the system and environment 
were encompassed in the initial state pse- We now take the opposite approach and now 
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start with a state ps of system S only. We then wish to embed this state into a state on 
the combined system. This is achieved by an assignment map: a map A from a density 
matrix of system S to a density matrix of the combined system SE given by 

A : p s i-> ,40s) = (4.2) 

where we still require that Tr^ps^) = p^. Any correlations between the system and 
environment are now a accounted for in the assignment map A. 

When using assignment maps in our description, any trace preserving quantum oper- 
ation £ on system A can be expressed as the composition of an assignment map A and a 
CP map T by 

£(p) = Tr E {ToA{p)). 
T is simply the unitary evolution of the composite system T (A(p)) = U (A(p)) U\ This 



is almost the same as Eqn. (3.10 ) used to describe the evolution in our system-environment 
model for CP quantum operations, however the assignment map means the state of the 
composite system is not restricted to being a product state. In this case £ is a CP map 
if and only if the assignment map A is CP |22j . 

The problem with assignment maps is that their definition is ambiguous, there is no 
restriction on how we embed our system into the combined system. This has been a cause 
of controversy in relation to using non-CP maps to describe physical processes [TH EE] • 
We need a method for characterizing which non-CP maps are physically relevant. 

A definition for a non-CP map £ being physically accessible [5] is that there must 
exist an assignment map A such that £{p) = Ti E [UA(p)U% where A need only act on 
a finite-volume subset of the state space of S. In addition a non-CP map only has a 
physical meaning when acting on its domain of positivity, where £(p) > 0. Physically 
relevant maps should also be able to be identified by quantum process tomography, which 
is a technique for characterizing quantum operations that we introduce in Chapter [5] 

Terno et al. [5] define two main classes of assignment maps which are physically acces- 
sible: Linear, and non-linear. These give rise to linear and non-linear system-environment 
correlations respectively. The linear case is the simplest scenario where the state of the 
environment is independent of the state of the system. For example A(p) = 7( p ) where 
Tts[7( p ))] = p but Trs[7( p )] = r Vp This type of assignment map will always lead to CP 
evolution as the state of the environment is independent of the state of the system [5]. 

In the non-linear case the state of the environment may be a function of the input 
state p. Not much is known about the non-linear case other than such assignment maps 
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often lead to non-CP dynamics [5]. We will show this with example in Section 6.2, when 
what is effectively a bilinear assignment results in a non-CP quantum operation. 



4.3 Quantum Discord 



Our example in Section 4.1.1 showed that even separable, or classically correlated, initial 
states can give rise to non-CP evolution. This suggests that taking classical and quan- 
tum correlations to be synonymous with separable and entangled states respectively is 
inadequate. A new method for defining quantum and classical correlations was proposed 
by Ollivier and Zurek [7], and independently by Henderson and Vedal [8]. This method, 
which we will now introduce, is known as quantum discord and uses two inequivalent 
quantum versions of the classically equivalent formulas for the mutual information that 
quantifies the correlations between two systems. We begin by briefly reviewing some 
classical information theory. 

4.3.1 Classical Information Theory 



In classical information theory [23J the Shannon entropy is a measure of the ignorance, 
or missing information, about a random variable A. It is given by 

H(A) = -J2p(A = o) \ogp(A = a), 

a 

where p(A = a) is the probability of A taking the value a. If we now consider two random 
variables A and B, the conditional entropy of A given B is 

H{A\B) = ^P(B = b)H(A\B = b) 
b 

= ^P{B = b) ^P{A = a\B = b)\ogP(A = a\B = b), 

b a 

where P(A = a\B = b) is the probability of A having value a given that we know B = b. 
The correlations between A and B are then measured by the mutual information 

J{A:B) = H{A)-H{A\B). (4.3) 



In all expressions the probabilities are derived from the joint probability distribution 
P(A,B). 
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It is possible to formulate an equivalent expression for Eqn. (4.3). By using Bayes 
theorem one can rewrite the conditional entropy as H(A\B) = H(A,B) — H(B) [23] . 
This gives us a second expression for the classical mutual information 

I(A:B) = H(A) + H(B)-H(A,B). (4.4) 

Our next step is to generalize the expressions J(A : B) and X(A : B) to quantum 
systems. 

4.3.2 Quantum Expressions for Mutual Information 

To generalize mutual information to quantum systems p], the random variables now 
represent the states of quantum systems A and B. The equivalent of the joint probability 
distribution is the density matrix of the state of the combined system, pab- The reduced 
density matrices pa = ^b(pab) and ps = Tta(pab) for systems A and B respectively 
are the equivalent expressions for the probability distributions of the individual systems. 

Our measure of ignorance about the state p of a system is then given by the von 
Neumann entropy 

S(p) = -Tr(plogp) (4.5) 
= ^ -A, log A, (4.6) 

i 

where {Aj} are the eigenvalues of p, and we define OlogO = if zero is an eigenvalue of p. 
The logarithms are taken to base 2. 



A quantum expression for Eqn. (4.4) is then obtained by replacing the Shannon en- 
tropies with their equivalent von- Neumann expressions. So the quantum mutual informa- 
tion is 

1a-.b{pab) = S(p A ) + S(p B ) - S(pab)- (4.7) 



In the case of Eqn. (4.3) however, the generalization is not trivial. In the case of 
conditional entropy, to express the state of system A given that we know the state of 
system B requires us to perform a set of measurements on system B. This is done by a 
complete set of 1-dimensional orthogonal projectors {ITf } acting on the system B. If the 
outcome corresponding to a measurement of Ilf is detected, the state of the joint system 
will be given by 

Pm? = -nf PAB iif , (4.8) 

Pi 

where pi = Tr(LTf pab) is the probability of the i th outcome. 
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The generalized expression for the conditional entropy is then obtained by averaging 
the entropy over the post-measurement states. That is 

S(p AB \{Tl?}) = J2p l S(PA\nz)- (4-9) 

i 

Our second quantum expression for mutual information gives a measure of the in- 
formation about system A which can be obtained by measuring system B. It is given 
by 

Ja:b(pab) {ii?} = S(pa) - S( PAB \{Uf }), (4.10) 
and depends not only on the joint system state pab, but also on the measurement set 

{nf}. 

The difference between the quantum generalizations of mutual information gives us a 
a measure of the non-classical correlations of a joint state. This difference depends on the 
choice of projectors used in J . In general we are interested in the minimum value of such 
a quantity, so we define the quantum discord to be minimum difference over all possible 
choices of projectors, 

T>a-.b(pab) = min [1a-.b(pab) - Ja-.b(pab)} ■ (4.11) 

S(pb) - S{p AB ) + ^PiS{p m B) 

i 

Some properties of quantum discord are 

1. Quantum discord is a non-negative quantity, i.e. T>a-.b(pab) > 0. 

2. For a state pab, ^a-.b(pab) = if and only if there exists a set of orthogonal 
projectors {nf } acting on system B such that 

Pab = U i P abU ? = J2 t >® n - ( 412 ) 

i i 

where n = Tr B [nf p^Ilf ] . 

Proofs of these results can be found in [7]. 

Now that we are done with the mathematical derivation of quantum discord we move 
on to discuss its physical significance. 



mm 

{nf} 
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4.3.3 Quantum Discord in Open Quantum Systems 

Quantum discord provides us with a quantitative measure of the non-classical correlations 
in a bipartite system. In a bipartite system, if we can obtain information from one system 
by performing measurements on the other then the systems are correlated. The extreme 
case of this is quantum entanglement. However, we can have non-entangled states which 
still exhibit this property characterized by a discord greater than zero. 
For example, consider a bipartite system of two qubits in the state 

P =\{Ph ® Pd + Pv ® Pa + Pd <%> Pv + Pa® Ph) , 

where pu = \H)(H\ , py = \V)(V\ Pd = \D){D\ and pa = |^4)(^4| • This state is separable, 
and the reduced density matrices of both systems are maximally mixed. We consider 
orthogonal projectors of the form ILf — I® where 

|Vo> = cos(0)|O) W^sm(0)|l>, 
= sin(0)|O) -e^cos(0)|l). 



A plot of Za-.b — Ja-.b of p as a function of 6 and <f is shown in Fig. (4.2(a)). One 
choice of projectors for which the quantity is minimum is 9 = f = 0. In this case we can 
calculate the discord to be 

V A .. B (p) = -\og (-\ ^ 0.33128. 



4 ° 4 



In the extreme case of a maximally entangled state, for example p^ = \4>){4>\ where 
= 4g(|00) + 1 11)), we have that T>a-.b(P4>) = 1 for all values of 6 and <f. This is shown 



in Fig. (4.2(b)) 



Ollivier & Zurek argue that quantum discord should be used to define what we 
mean by classical and quantum correlations. The term classical correlations should apply 
to systems with zero quantum discord, and quantum correlations to systems with QD > 0. 
In their paper they reverse the order of the systems in their definition of quantum discord 
when consider a bipartite open quantum system AB consisting of a principle system A 
interacting with an environment B. This means the discord is given by 

D B :a(pab) = min (1 b -.a{pab) - Jb-.a(pab)) , (4.13) 
{nf} 

where the set of projectors {Hf} being minimized over act on system A. In addition we 
now have Vb-a(pab) = if and only if there exists a set of orthogonal projectors {Hf} 
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1.5^0 

(a) Separable state p 



(b) Maximally entangled state \4>){<t>\ 



Figure 4.2: Za-.b — Ja-.b a function of projectors for the separable state p 
\(ph <8> Pd + Pv ® PA + Pd ® Pv + Pa® Ph) and a maximally entangled Bell state \<j>) 
^(100) + |11)). 



acting on system A such that 

P AB = nfPABllt = U i ® ( 4 ' 14 ) 

We shall differentiate between the two expressions with the subscripts A : B and I? : A. 
Their argument for this is that in general one does not have access to the environment, 
and thus measurements should be performed over the system A. 

The main result of [6] is the following theorem. 

Theorem 4.1 The reduced evolution of any classically correlated bipartite system 
(Db-.a(pab) = 0) is always completely positive. 

The original proof for this is given in [B], however the notation used is not particularly 
clear. We produce our own proof for this theorem in Appendix |A.2.1 in what we believe 
to be clearer notation. 

The presence of quantum correlations allows for the possibility of non-CP quantum 
operations. However, while the presence of quantum correlations is found to be necessary 
for non-CP dynamics, it is not sufficient. While this provides a convenient method for 
determining situations where CP evolution is guaranteed, and extends the known class 
of systems which will always have CP evolution from simply separable to those with 
classical correlations, it does not help us deal with systems which may potential have 
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non-CP evolution. We have found some interesting issues which arise from the fact that 
quantum discord is not symmetric. 



4.3.4 Asymmetry of Quantum Discord 

The definition of quantum discord in the literature is ambiguous. In Ollivier and Zurek's 
original derivation [7] we are simply considering two systems A and B, the choice of 
orthogonal projectors used in generalizing the condition entropy to quantum systems has 



the measurements carried out over system B, which gives us the definition of Eqn. (4.11). 



This approach has been used in several other following studies j2H [25] - Alternatively we 
could take the approach used by Sudarshan et al. [6] reversing the order of the systems 



in the definition giving us an expression for discord as defined by Eqn. (4.13) 



The issue here is that quantum discord is not a symmetric quantity, in general T>ab{p) 7^ 
T^b-.a{p)- More importantly, a discord of zero in one direction does not imply zero discord 
in the reverse. We will show this by a counter example. 

Suppose T>a-.b(p AB ) = 0. This implies that there exists a set of orthogonal projectors 
{nf } on B satisfying Eqn. (4.12). If this also implied that T> b-a{p ) — 0, then there must 



also exist a set of orthogonal projects {Hf} acting on system A which satisfy Eqn (4.14). 
Consider the case where p AB = | (pn ® Ph + Pd ® Pv)- It is easy to see the orthogonal 
projectors ILf = p#,nf = p v satisfy Eqn. (4.12) and hence T>a : b(p AB ) = 0. 



Now suppose there exists Iky , Ilf satisfying Eqn (4.14). For this to be true we require 



that 



U A p H U A + U A p H U A 
U A p D U A + U A p D U A 



PH, 
PD- 



(4.15) 
(4.16) 



However pu and p^ are not orthogonal, so such a set of projectors cannot exist, hence 
T^b-.a{p AB ) 7^ 0. In fact a numerical computation gives the value to be T>b-.a(p AB ) = 



0.2018. Fig. (4.3) shows how T>a-b and T>b.a vary as function of the projectors used. Here 



we can also see that the distributions are different. 

This presents us with a particular problem in the physical interpretation of quantum 
discord. If we are to call two systems classically correlated when there is zero discord 
between them, which definition do we use? One motivation, which is implicit in the 
interpretation of Sudarshan et al. [6] , is that in the case of an open system A interacting 
with an environment B, we can not in general make any measurements of the environment. 
Thus the definition with projective measurements acting on the system gives us a measure 
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of how much information we can deduce about the environment by only measuring the 
system. 

An alternative to this approach is to ask the question how much information can 
the environment know about the system? In this case it is the reverse T>a-.b which is 
important, as it reveals how much information about A can be obtained by measurement 
on its environment. 

Each approach is supported by a valid argument and thus we are left with a unique 
problem. In the case of our example, according to T>a-.b = we have a classically correlated 
state, yet T>b.a 7^ implies that the state in fact has quantum correlations. Our state 
should be either one or the other, unless correlations in quantum mechanics are directional. 



An important implication of this is in generalizing Thm. (4.1). 
We propose the following conjecture: 

Conjecture 4.2 If an initial state p of a bipartite system AB is classically correlated, as 
defined by T>a-.b(p) = 0, yet also has T>b-.a(p) 7^ (so the correlations are not symmetric), 
its evolution need not be completely positive. 

Our basis for this conjecture is as follows: We try and follow a similar argument to 



the proof for Thm. (4.1) from Appendix A. 2.1 but using the condition from Eqn. 4.12 
for Va-Ap) = 0, 
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Since the projectors IT are orthogonal pure states we can write them as IT = 
Evolution of the joint system by an arbitrary unitary operator U then leads us to 

£{p) = Uprf = J2p* ( ^(nfUlmfSmi J n (j2S mj (m\ B E7+ \n) B 

i,j,n \ m / \ m 

Now, for each there exists a set of projectors Xjj such that t { = ^\ ^i,j T i^i,j — 
5^- qijXij. In general each X^j are not orthogonal and do not commute for different i. 
This means we use the properties of these orthogonal projectors to remove the dependence 
of the Kraus operators on i as we did in the previous case. Hence it appears that in general 
the Kraus operators will be different for different states of the reduced system, suggesting 



that the evolution in general is not completely positive. Thus Thm. (4.1 ) is only true for 
T^b-.a{p) = 0, it does not appear to generalize to the case of T>A:b(p) = 0. 

This indicates that the definition of classical and quantum correlations based on quan- 
tum discord is insufficient, as it is not symmetric with respect to the two systems A and 
B. A possible avenue for future research into this result could be to consider what would 
happen to the definition of quantum discord if we allowed for more arbitrary measure- 
ment sets, such as projective operator valued measurements (POVM), or to consider a 
symmetrized version of the quantum discord. 



5 

Quantum Process Tomography 



We have now introduced the mathematics needed for describing evolution of open quan- 
tum systems by completely-positive trace-preserving maps, and also some methods of 
characterizing correlations of a joint system state. It was shown how certain initial cor- 
relations of an open system with its environment could lead to non-completely positive 
evolution. In this chapter we will now introduce a technique for characterizing unknown 
quantum operations on open systems. This technique is called quantum process tomog- 
raphy and allows for the characterization of an unknown quantum process by measuring 
the output states for a relatively small known set of input states and using this data to 
reconstruct a process matrix for the evolution. 

Several different schemes have been proposed to perform this task. We will investigate 
two of them in this chapter. The techniques are standard quantum process tomography, 
and ancilla assisted process tomography. But first, before we can properly describe the 
reconstruction of an unknown quantum process, we must introduce a method for recon- 
structing unknown quantum states. 
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5.1 Quantum State Tomography 

If we have a quantum system in an unknown state we require a technique to completely 
determine the density matrix which represents that state. Such as technique is quantum 
state tomography. 

In state tomography an experimenter reconstructs the density matrix by inverting 
probability data obtained from measuring many identical copies of the unknown state 
with an appropriate set of measurements. Such a set of measurement operators, {Mi}, 
must have the property that any density matrix of the system can be expressed as a linear 
combination of these operators, i.e. p = ^ a^Mj. Such a set is called tomographically 
complete, and forms a basis for the density matrices of the system. 

In the case of a single qubit represented by the polarization state of a photon, a 
commonly used tomographically complete measurement set is {Mj = where i = 

H, V, D, R. These are the density matrices corresponding to the horizontal, vertical, 
diagonal, and right-circular polarization states of a photon respectively. Even though 
Mi = pi, as pure state density matrices are projective measurement operators, we will use 
M to distinguish between states and measurements. 



Recall from Section 2J3 that for an unknown state p, the probability of detecting the 
outcome corresponding to measurement Mj is given by Pi = Tr (pMj). With an appropriate 
basis {Di} we can reconstruct the density matrix p from these probabilities by [26J 

p = Y t p i D i . (5.1) 

i 

The basis {D,} used for reconstruction is called the dual basis for {Mi}. It is defined 
by the orthogonality relationship 

Tr (ujMy) = Sij. (5.2) 
Calculations for the dual basis for input basis {Mh, My, Md, Mr} can be found in Ap- 



pendix A.4 For systems of multiple qubits we can form the tensor products of the mea- 
surement set for a single qubit and do joint measurements. For example, if p is a 2-qubit 
state a tomographically complete measurement set will be {My} where My = Mj ® My 
Since measurement in quantum mechanics affects the state of the system, we require 
many identical copies of an unknown state in order to accurately determine the proba- 
bilities Pi. If, for each measurement operation Mj, we prepare N copies of a state p, the 
number of detections recorded, iVj, allows us to compute the probability coefficients in 
the standard way. 
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Since we are determining the relative probabilities by count statistics, statistical noise 
becomes a factor in reconstruction. This can lead to the reconstruction of unphysical 
states. In these cases optimization techniques are employed to optimize the reconstructed 



state to the closest physical one. This process is discussed in more detail in Section 5.3.1 



5.2 Quantum Process Tomography 

Now that we know how to characterize unknown quantum states we can move on to the 
characterization of quantum processes. First we shall describe the general motivation 
common to the various techniques of tomography. 



Recall from Section 3^ that the evolution of an open quantum system is commonly 
described by a quantum operation: a linear, trace-preserving, completely positive map £ 
acting on the state space of the system. There were several mathematical representations 



presented in Section 3.2 which allowed us to completely describe the evolution of an 
arbitrary state p by such a CPTP map. 

Now we consider the reverse situation. If we are presented with an unknown operation, 
but believe that it is CP, what technique might we employ to characterize it? One option 
is to send in known quantum states and measure their outputs. However, it is a practical 
impossibility to do this for all possible input states. A characterization technique which 
requires only a finite set of inputs is process tomography. 

In general we wish to find a mathematical representation for £ through experimental 
observations. The methods we will describe in the following sections are formulated in 
terms of the process matrix A5 which uniquely describes £. The two methods we will 
consider are called standard quantum process tomography (SQPT), and ancilla assisted 
process tomography (AAPT). We will begin with SQPT. 

5.2.1 Standard Quantum Process Tomography 

One of the earliest proposed methods for performing process tomography is what is now 
referred to as standard quantum process tomography (SQPT) [27] . The idea of this scheme 



is shown in Fig. (5.1). An experimenter prepares a basis of input states, subjects them 
to an unknown operation £ and determines the output states by state tomography. This 
information is then used to reconstruct a process matrix for £. 

We start by choosing a basis of density matrices of the principal system {pj}- This 
means that any density matrix p can be uniquely expressed as a linear combination 
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{M k } 



Figure 5.1: Schematic of standard quantum process tomography 

p = ^2jPjPj- For a rf-dimensional system, we require a set of d 2 linearly independent 
density matrices. 

The choice of an input basis is not unique, and is often determined by experimental 
convenience. For single qubits in linear quantum optics experiments it is common to 
choose the states corresponding to horizontal, vertical, diagonal, and right-circular polar- 
ized photons. I.e. the basis {\H) , \ V) , \D) , \R)}. The corresponding density matrices for 
these states are: 



Ph 
Pd 



(I 



0"3 



Pv 
PR 



(I - C7 3 ) 
(I+<7 2 ). 



Notice that we chose the same basis for input states as for the state tomography mea- 
surement operators. 

To construct the process matrix we use the approach of Kuah et al. [9]. However, their 
notation is confusing and hard to generalize. We will formulate an equivalent description 
inline with the mathematics introduced in Chapter [3] 

Let 8 be an unknown quantum operation acting on an open quantum system. We 
will denote our basis of input states for the system by {pj}. The output states will be 
denoted by {£(/??■)}. 

We can then define a process matrix for our quantum operation by 



(5.3) 



where Dj is the element of the dual basis corresponding to pj, and Dj designates the 



complex conjugate of the elements of Dj. See Appendix A. 5 for the proof that this 
expression is equivalent to our original definition of the process matrix. 

The output states are determined by performing state tomography as outlined in 
Section [5TT If we use the same basis for our input states and measurement set we can 



combine Eqns. (5.1) and (5.3) to give 



5.2 Quantum Process Tomography 



41 



A £ = J2MpiM J )(D*®D J ), (5.4) 

ij 

and if our input and measurement sets differ each will just have a different dual basis. So 
the total number of different input state and measurement combinations is {d 2 ) 2 . 

From here we can use the equivalences between different representations of CP maps 
described in Section [372 to convert to a Kraus representation, or superoperator if such a 
representation is preferred. 



5.2.2 Ancilla Assisted Process Tomography 



Another method of process tomography is ancilla assisted process tomography (AAPT) [28l 
29] . This method is mathematically equivalent to SQPT, and exploits the equivalence 
between quantum operations and quantum states of a larger system we introduced in 
Section 1372721 

Recall from the Jamiolkowski isomorphism that for a quantum operation £ acting on 
a (i-dimensional system S, the process matrix could be defined as 



A^ 
d 



(J®£)(|0)(0|) = P£, 



(5.5) 



where \<j>) = <8> \j))/yd is a maximally entangled state of the joint system AS 

where A is an ancilla of the same dimension as S. 

Hence if we introduce an ancilla A and prepare the maximally entangled state on AS, 
we subject only half of this joint state to the operation £ and determine the output state 



ps by state tomography. This is shown in Fig (5.2). We them obtain the process matrix 
simply by A £ = dp £ . 




DW) 



Figure 5.2: Schematic of ancilla assisted process tomography 
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AAPT has the advantage of only requiring the preparation of a single input state 
instead of the d 2 required for an d-dimension system in SQPT. However, since the input 
state must be a joint state of our principal system and an ancilla of equal dimension, the 
required number of measurements to perform state tomography on the output is increased. 
At the end the total number of measurements is the same as in the case of SQPT, {d 2 ) 2 , 
as we must measure on the larger joint state space. 



AAPT With Alternative Input States 



Eqn. (5.5) is only valid for the maximally entangled input state |0)(0|. Suppose we want 
to use another pure input state described by a state vector \M) = m ij K) ® \ Then 
using our relationship between vectors and matrices from Section |2.10[ we have 

\M) = VdM T ®I|0) = Vdl ® M \(f)) , 

where M = mat \M) = Y2ij m ij Ij'X^Ij anc ^ 10) 1^) l^fd- 

Thus the output state generated by X® £ and the pure input state pu = \M)(M\ is 

{X®S){p M ) = {T®£){\M){M\) 

= (VdM T ®l) [{I®£) |0) {VdM* (8)1) 

= (M T ®l)dp £ (M* ®l) 

= (M T ®1)A £ (M* (g)I). 

Hence we can recover the process matrix by 

A £ = (M T )- X ® I [(£ ® T)p M \ (M*)- 1 ® I (5.6) 

This is possible if and only if M is invertible, ie. if and only if rank(M) = d. In terms of 
the vectorized matrix this is equivalent to requiring Sch \M) = d, which means \M) must 
be entangled. 



AAPT With Separable Input States 

The method of AAPT just described is also referred to as entanglement assisted process 
tomography due to the input being an entangled state. What happens if we wish to use a 
separable input state p? By the spectral decomposition, we can express this as the sum 
of pure states p = J2k \Mk){Mk\ where \Mk) = J2ij(^k)ij \j) <%> The output state of 
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p generated by £ is then 

{1®£){ P ) = ^(Mj®I)A £ (M*®I) 

k 

= (M®l)(A e ), (5.7) 



where M(p) = J2 k M kP M k, so the 

process matrix is recoverable through Eqn. (|5. Th if 
and only if the map M. is invertible. An input state satisfying this is called faithful. 

While we need not use an entangled input state for AAPT it has been suggested in 
literature that these states are the most efficient [2H]- Here efficiency refers to requiring a 
smaller ensemble of copies of our input state to accurately determine the output though 
state tomography. This is due to statistical errors which result in quantum measurement, 



the nature of which are discussed in Section 15.11 

A measure of the faithfulness of an input state pab with maximal Schmidt number is 
given by [30] 

F{ PAB ) = Tr(p\) = ^ A- 



\ 2 
s ■ 



where p a = Ttb[pab} and Aj are the eigenvalues of pA- This coincides with the definition of 
purity in Section |2.2 The only pure states with maximal Schmidt number are entangled. 



5.3 Non- Completely Positive Maps in Process To- 
mography 



We mentioned in Section |5.1| that non-physical output states often occur as a result of 
statistical noise in process state tomography. The effect of these reconstruction errors can 
lead to non-CP maps being observed in process tomography. The standard formalism for 
process tomography is to assume that evolution must be completely positive. Then any 
deviations from complete positivity in our output state must be the result of statistical 
and/or experimental errors. In this case a numerical optimization of the output state is 
done to enforce the requirement of a CP process map [3TJ [321 [33]. A common imple- 
mentation of this idea is maximum likelihood process tomography, which we will discuss 
below. 

However, there is a problem with this approach. As we have previously mentioned, and 
indeed showed by example in Section 4.1.1] non-CP maps can arise from the presence of 



initial correlations between a system and its environment. If we can have situations where 
a non-CP process map should be observed through tomography, employing optimization 
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techniques based on the assumption of complete positivity will give us an incorrect process 
matrix. This leaves us with a situation where we must distinguish between a non-CP 
process map arising due to statistical noise, and one arising due to initial correlations. 
This question will be the focus of Chapters [6] and [7| 



5.3.1 Maximum Likelihood Process Tomography 

Before we move onto a discussion of state preparation in process tomography we will 
briefly describe the basic theory of maximum likelihood process tomography. This will 
be important when we investigate the statistical distribution of reconstructed process 
matrices in Chapter [7j 

In SQPT, for a given process map and input states {pj}, the expected number of 
counts for M measurements of an operator M m is 



NTi A [(pj <g> M m )A e ] . (5i 



For optical experiments such as we will consider in Chapter |8j the photon sources are 
typically spontaneous. This means they will produce photonic qubits at random intervals, 
though at an overall average rate. The count statistics in this case will be Poissonian [31]. 
If we assume that the length of the experiment is sufficiently long, we can approximate the 
Poissonian distribution by a Gaussian distribution, with mean and variance both equal 
to rij m [35]. Under these assumptions, the probability of detecting a sequence of counts 
{ n )m\ is g iven by 



P({n e jm }\A £ ) cxHexv 



,/»' 



n e jm - n jm (A e )) 



2n jm (A £ ) 



(5.9) 



where n jm is the expected number of counts given by Eqn. (5.8), and P and rijk are 
functions of Ag. 



Eqn. (5.9) is called the likelihood function, and we wish to maximize this function 
over the variable A& to find the process matrix most likely to give the detected sequence 
of counts. This is maximum likelihood estimation. However, it is generally easier to 
minimize the negative of the log-likelihood function, so the actual optimization task is 

K m - n jm (A £ )] 2 



minimize 



4^ n jm (A e ) 
with constraints A^ > and Tr^A^) = I. 
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It should be noted that optimizing this log-likelihood function is equivalent to a 
weighted least squares problem, however the weights depend on the process matrix. In 
addition the requirement that the process matrix correspond to a physical state give us 
our constraints. 

It is possible to write this seemingly nonlinear optimization routine in the form of con- 
vex optimization called a semi- definite program [36], for which there are efficient numerical 
tools for solving [3"T] . 
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State Preparation In Process Tomography 



To perform quantum process tomography, as outlined in Chapter [5j the use of specific 
input states was required. However, we made no assumptions as to how these states were 
prepared. In this chapter we investigate methods for preparing an arbitrary state of an 
open quantum system into the required input states for a tomography experiment. It is a 
particularly important issue if a system and environment are initially correlated since the 
evolution of the system may not always be completely positive. We investigate how state 
preparation can influence the complete positivity of system's evolution. This approach 
is motivated by the work of Kuah et al. [9] which investigated two types of preparation 
procedures in SQPT, but both within the context of CP dynamics. We use their paper 
as a starting point to investigate the situations when non-CP maps can arise. 



6.1 State Preparation in SQPT 

As mentioned in Chapter [5j an unknown quantum operation £ can be completely char- 
acterized by the use of SQPT. To perform SQPT we require the preparation of a select 



set of input states. The idea of this procedure is shown in Fig. (6.1). If we have an open 
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system S with environment E, the initial state of the joint system, 70, may be correlated. 
State preparation is a procedure applied only to system S to prepare 70 into the required 
tomographic inputs {pj}- Following preparation, SQPT is performed as previously de- 
scribed. We will now review methods of state preparation proposed in the literature, 
and also propose our own scheme, and investigate how they fare in the presence of initial 
correlations. 



7o< 




E ' 



Figure 6.1: State preparation in SQPT 



6.1.1 Stochastic State Preparation 

First we consider the so called stochastic state preparation j9]. The general idea of this 
procedure is to remove correlations by a process which sends the state of the principal 
system S to a fixed state regardless of the initial state of the joint system SE. For example, 
cooling a quantum system will bring it into its ground state. From here stochastic unitary 
operations are used to transform the produced state into the set of input states needed 
for SQPT. 

Mathematically we describe this procedure in terms of a preparation map. The first 
step of this procedure is represented by a map 5 acting on density matrices of system S 
which sends all states to a fixed pure state \ip){ip\. So 

H(p) = \4>)(4>\ Vp of system S. (6.1) 

Applying this map to states of the joint system gives 

(H® J)(7) = IV'XVI <8> t 

where r = Tr^ [(H (8)2") (7)] is the post-preparation state of the environment. 



6.1 State Preparation in SQPT 
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With joint system in the product state <8>r, the next step is to prepare the basis 

of input states pj needed to for SQPT. This is done with a set {Sj} of so called stochastic 
maps which satisfy Sj(\ip)(ip\) = pj. In fact these are simply a set of unitary operators 
chosen to give the required states pj when applied to the fixed output state of H(p). 

Hence the stochastic preparation procedure is given by the collection of maps {Pf°}, 
where 

Vf° = ( l S j oH)®X (6.2) 
Vf°(a) = pj®r Vp of the joint system SE. (6.3) 

If our system is initially in the (possibly correlated) state 70, and undergoes unitary 
evolution via a quantum gate U, SQPT performed on the output states 

S( Pj ) = Tr B [UVf ( 7o )C/t] 
= Tr B [U(j> j ®T)U*]. 

We note that the authors [9] claim the state r is a function of the preparation map 
H alone, hence it will be fixed for any input 7. However, if the map H is itself trace 
preserving this is incorrect. In this case r is independent of the map S and only depends 



on the initial state 7. We prove this in Appendix A. 6.1 Since the state of the environment 



t is fixed for all pj, by Thm. (3.3) the resulting process matrix from SQPT will always 
be CP. 

6.1.2 State Preparation by Measurement and Rotations 

We now propose a preparation procedure similar to the stochastic case which could be 
implemented in a linear optics SQPT experiment. In [9] the authors consider using a set 
of projective measurements as a preparation procedure. We will discuss the procedure 
using only measurements momentarily, but first we will introduce our own procedure 
based on this idea. We consider replacing the map H from the stochastic procedure with 



a single projective measurement n. Recall from Section 2.5 that in quantum mechanics 
measurement is described by a set of measurement operations Mj. If we measure a state 
p and detect outcome j, then it is assumed the post-measurement state of our system is 
Pj = ^-MjpMj, where p = Tr(pn) is the probability of detection. 

We consider the case of measuring with only a single projective measurement operator 
n. For our example we will use Pin = \H)(H\ which projects onto the horizontal polar- 
ization state ph- This could be implemented in linear optics for example by a polarization 
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filter set to only let through horizontally polarized light. If a beam of randomly polarized 
light is incident on the filter, only photons in a polarization state corresponding to the 
filter will pass through. The sacrifice in this is that we are performing a trace decreasing 
operation. We will loose qubits orthogonal to the projection. This means we may have 
to judiciously select our initial projector based on the initial state of the joint system so 
that we are not loosing a large portion of our qubits. 

In any case, once we have performed the projective measurement our stochastic maps 
are then the set of unitary rotation matrices {Rj} such that RjPhR] = Pj where pj is one 
of our SQPT input states (if, V, D, R). In optics such rotations to the polarization state 
of a photon can be easily implemented by wave plates. 

Our preparation maps are then given by 



V 



f°(7o) = ^{R^u) A l^uR\) A (6.4) 
= Pj ® Tff (6.5) 

where {RjU.H) A = (R^h) <S> I is an operator acting only on the first system, r# = 
Tr A (RjU H ) A l0 (U H R]) A /r H is the post-preparation state of the environment, and 
Tji = Tr [n^7o] is the probability of detecting the outcome corresponding Uh- Here 
we see that in general the post-preparation state of the environment is dependent on both 
the initial state 7 and the measurement operator Ii H . 



6.1.3 State Preparation by Projective Measurements 

In our example of a stochastic preparation procedure we introduced the idea of using 
a projective measurement followed by unitary rotations to prepare the state of our sys- 
tem. We now propose the idea of using only measurements to prepare our input states 
for SQPT. This method of preparation would be convenient in linear optics where such 
measurements can be implemented with polarizing elements. 

If {pj} is our required input set for SQPT, we chose a set of projective measurements 
{Hj} where n,, = pj. Hence the preparation procedure for projective measurements is 
given by the collection of maps {p™ eas | where 

Vr as (l) = ^-nf-yUf (6.6) 
1 3 

= Pj®Tj. (6.7) 

Here Tj = Tr A [n^nj 1 ] /Tj is the post-preparation state of the environment, and Tj = 
Tr [11^7] is the probability of detecting the outcome corresponding Tlj. 



6.1 State Preparation in SQPT 
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Since we are using measurements alone for preparation we have the possibility of 
failure in our preparation. If the initial state is orthogonal to one of our projectors IT,-, 
then we can never prepare the required state pj. If the initial state is instead only close to 
orthogonal, we will have a very small probability of producing the required state. Since 
we must prepare many copies of the state for performing state tomography this greatly 
reduces the efficiency of the experiment. 

If our system is initially in a state 70, and undergoes unitary evolution via a quantum 
gate U, the output states for our set of tomography input states are given by 

£( Pj ) = Tr B (UP™ as ( l0 )tf) (6.8) 
= Tr B (UpjQTjtf). (6.9) 

Here we can see that in general the state of the environment, Tj, depends on the state 
of the system, pj. Since the state of the environment is not fixed, the process matrix 
determined by SQPT will not necessarily be CP. This situation can occur when we use a 



correlated initial state 70. We will show this with a simple example in Section 6.2 



6.1.4 Non-CP Maps due to Preparation by Measurement 

We now ask whether we can find some bounds on the type of correlations which may give 
rise to non-CP process matrix. Our initial idea was to use the concept of quantum discord 



which we introduced in Section 4.3 Recall from Thm. (4.1) that the evolution of any 
initial state 7 satisfying Ps:s(7) = will always be completely positive. We find that if 
one uses projective measurements alone to prepare the input states for SQPT, then this 
theorem no longer holds. An initial state 7 satisfying T>e-.s{i) — can lead to non-CP 
evolution in SQPT. 

To see why this occurs consider an initial state 7 = | (ps <8> pu + Pv ® Pd) which, by 
Eqn. (4.14), satisfies 'Db:a{i) = 0. Now if we use projective measurements to prepare a 
set of input states {pi = where % = H,V,D,R, the post-preparation states of the 

environment are 



t h = Ph 

T V = PD 

T D = ^(pH + pD) 

tr = ^(Ph + Pd)- 
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Even though the preparation procedure puts SE into a product state, since the state 
of the environment is not fixed, the evolution need not be CP. 

The situation for ensuring that the evolution of a state 7 will always be CP when using 
SQPT and preparation by projective measurements is more restrictive than normal. Due 
to the nature of the preparation procedure, we not only require that the initial state be 
simply separable, 7 = ps <E> t e , we also require that the reduced state of the system, p s , 
not be orthogonal to any of our required input states. 



6.1.5 Non-Linearities from Preparation by Measurement 

Our next issue with preparation by measurement is to do with the non-linearity of the 
procedure. In [9] the authors raise the issue of non-linearity with this preparation proce- 
dure. While the process matrix determined by SQPT from this procedure will be linear, 
the preparation procedure itself is not. 

For example, if we wished to prepare the system into a state p = \{pi+ P2), where 
and pj are elements of our input set, then we have 



spmeas/^ 



Til + n 2 ) A 7 (n 1 + u 2 ) A 
4Tr [i(n 1 + n 2 )^ 7 ] 



2(ri + r 2 



So the preparation procedure is in fact bilinear. A method of dealing with this situation 
by describing the evolution with a bilinear process matrix is proposed [9] . The method of 
using a bilinear process matrix is significantly more complicated than the linear case, and 



we discuss this approach in section (6.4). Another problem with this assumption is that 
using projective measurements, we can only ever prepare our system into a pure state, as 
mixed states do not satisfy p 2 = p. 

We propose an alternative protocol of preparation by measurement to ensure a lin- 
ear procedure. Our set of preparation maps | , pj neas } correspond to a basis for density 
matrices. If we wish to prepare a state p = J2iPiPi> where are the basis states, we 
simply perform the stochastic preparation procedure V meas [^i) = 'Yl li PiP^ neas { r )). This is 
done by randomly performing each of the basis preparation procedures with probability 
Pi. This is analogous to the preparation of mixed states by randomly preparing one of 
several different pure states weighted by their respective probabilities. 



6.2 Examples of Preparation Procedures in SQPT 
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It should be noted that this in general will leave the joint system in a state V meas {^j) = 
YliPiPi ® T ii which will still in general give rise to non-CP dynamics. 



6.2 Examples of Preparation Procedures in SQPT 



We will now consider an explicit example of performing SQPT with the two preparation 
procedures mentioned in Sections 6.1.2| and |6. 1 .3 We will consider a system S and its 



environment E, both qubits, with an interaction given by a controlled-NOT gate. This 
gate performs a bit flip operation (|0) <-> |1)) of the target qubit if the control qubit is in 
a state |1). If the control qubit is in state |0) the gate does nothing. We will let S be the 
target qubit, and E the control. The unitary matrix for this gate is given by 



U, 



( 


1 








o\ 













1 










1 










1 





0/ 



CNOT 



Suppose that the initial state of joint system SE is the maximally entangled state 
P<t> = 10) (01 where |0) = ( 1 00) + 1 11))/ v2- Our tomographic input states will given by pi 
where i = H,V,D,R. 

We now consider the resulting process matrix when we preform SQPT using the two 



preparation methods described in Sections 6.1.2 and 6.1.3 



• Method A: We use a single projective measurement Uh to prepare system A to 
the state pn, and then use unitary operations Ri to rotate pu to the state pi. 

• Method B: We use a set of projective measurement IT to prepare each of the initial 

states pi. 

The process matrix reconstructed though SQPT are 

10 1 





10 1 



A, 



and Kb 





2 





-1 - i 


1 


1 








1 


l + i 


2 




1 


2 







1 


1 -i 









(6.10) 



for method A, and method B respectively. 

has eigenvalues (2, 0, 0, 0), and hence Sa the process reconstructed from preparation 
by a measurement followed by rotations, is completely positive. on the other hand 
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has eigenvalues (1 + ^,—^,^,1 — ^) and hence £b the process reconstructed from 
preparation by measurements alone, is not completely positive. 

In the first case the evolution is not very interesting, it is simply the identity map 
Ea = I-, this is because the state of the environment after the preparation procedure is 
|0)(0|. 

In the second case we have a genuine non-CP map. To see why this happens let us 
consider the state of the environment after each preparation procedure. We have that 







So in each case the state of the environment is different. For example, in the case of 
projecting onto pu the state of the control qubit th means the evolution of qubit A will 
be the identity, U = I, while in the case of projecting onto pv our gate becomes a bit flip, 
U = X. Even though our initial state is simply separable, the state of its environment is 
dependent on the state of the system due to the nature of our preparation procedure. 

This is an interesting result, as conceptually the only difference between the two 
preparation procedures was an extra rotation stage that was introduced in Method A. 
This means that in a linear optics tomography experiment, having a polarizer to prepare 
states can result in a dramatically different result than a polarizer with a wave plate placed 
after it. Even though these elements act only on system S, the initial correlations between 
S and E give rise to a CP result in one case, yet a non-CP result in the other. This is 
explained in more detail in Chapter [8] where we propose an experimental implementation 
to illustrate this observation. 

While we used a maximally entangled input state, this result holds even for separable 
inputs. For example if we use the input state p = \ (p# ® Pa + Pd ® Pv), then we get 
a process matrix with eigenvalues {1.642,0.507,-0.253,0.105}, so the evolution is still 
non-CP. In this case the quantum discord for p is Va-.b(p) = T^b-.a = 0.1443. 



6.3 State Preparation in AAPT 

We now discuss state preparation in AAPT. This situation is much simpler than in SQPT 
as we only need a single input state. The difference is now we are dealing with three 



systems, an ancilla A, principal system S and environment E as shown in Fig. (6.2). We 
will assume that the joint system ASE is initially in a tripartite state jase- 



6.3 State Preparation in A APT 
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To prepare our input state for AAPT we perform a preparation procedure on the 



ancilla and principle system analogous to the 5 map from Section |6.1.1| This procedure 
prepares AS* in the maximally entangled state p^ s = \(j>)(4>\, where \<j>) = ^ \i) ® \i) /Vd, 
and is given by a map 



2(7as) = 10) (01 V7as of the joint system AS. 
Hence the preparation procedure on our whole system is given by 



(6.11) 



ASE) 



>Z){1ase) 



T E , 



where te = Tyab [(H <%>Z)jase\ is the post preparation state of the environment. This 
preparation procedure could be implemented by a Bell-state measurement. This is a joint 
measurement over two systems, which detects the maximally entangled states [30J . 



lASE (- 




\ E 



Figure 6.2: State preparation in SQPT 
Since the interactions between our system and environment are described by a unitary 



operation U, as shown in Fig. (6.2), the process map for this interaction, as described by 
the output state ps from AAPT is then given by 

ps = Tte[{1®U)V{ 1A se){1®U ] )\ 
= Tr E [{I®U){pf ®T E )(l®tf)\ 

Since we only need prepare this one input state to perform AAPT (in practice we could 



use another initial state so long as it satisfied the requirements in Section 5.2.2 the state 
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of the environment is fixed. Thus performing AAPT will be equivalent to performing 



SQPT with a fixed environment te and so, by Thm. |3.3[ the process matrix for this 
evolution will always be CP. 



6.4 Bilinear Quantum Process Tomography 

Even when we have a situation where preparation by projective measurements gives a 
CP results (for example when either direction discord is zero), if we treat the preparation 
procedure as an assignment map it will still be be bilinear. 

A method of performing bilinear process tomography was proposed in We will 
briefly present this method, and as done previously we will also convert their formulas 
into more familiar notation. 



We consider the same situation as section (6.1.3 ), where we have a system A interacting 



with environment B, both of dimension d. The initial state and evolution of the joint 
system being given by 70, and U respectively. If we project the initial state of the system 
onto p n for our preparation procedure, then the output is given by 

£{pV>) = ^ Tr B [U{pM ® J) 7 o(p (n) ® I)U^] . 

The idea proposed in [U] is to express this in terms of matrix elements, and then 
rearrange the elements so that 

£(p (n) )r, S = U re,r'a P^r" 10r"a,s"P P^l U l'f3,se 

r',r" s',s" a,fj,e 

= p(n) P r"l> ( U re,r>a 70r"a, S "/3 ^J'ftse J P s"l> 

r',r" s',s" \a,/3,e / 

- _L_ p(")* p( n ) 

~~ p(n) / t / t r r"r' Jvl r"r';s"s' r s"s'i 



r\r" s\s" 



where we have defined a new matrix A4 by 

M.r"r';s"s' = 5^ ^re,r'alOr"a,s" pU\, ^ Be . (6-12) 

The matrix M. to called the bilinear process matrix, and we can see it depends on both 
the initial joint state 70 and the evolution U. 



6.4 Bilinear Quantum Process Tomography 
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This representation of M. in is not particularly useful so we propose a more compu- 
tationally convenient one. Our idea is to break the d 2 x d 2 matrices U and 70 into d x d 
block matrices, where each block element is a d x d matrix. I.e 



/ 



U 



P\ 



12 



[U] 



Id 



\ 



\ P\ 



dl 



PU J 



13 Jm,n 



u, 



where ([£/] 
noting that [U% = [U]* 



(i—l)d+m,(j—l 



)d+n- Similarly we form the block matrices for 70 and W 



In this notation we can calculate each entry of M. by 

M { ;f)^ s , = TT([U} T Alo]riiAU}l, 



(6.13) 



M. itself is a block matrix, where (r'V; s'V) denotes the block elements of M, and (r, s) 
are the comments of a given block. Evolution is then given by 

1 



Qn 



(p m \ M \p n ) 



(6.14) 



where we have vectorized p, and each element of p acts on a block matrix element of A4. 
Verification of these expressions can be found in Appendix A. 6.2 



It was shown in [9] that M. is hermitian. However it was mistakenly asserted that 
Tr(A^) = 1. We found that Tr(A^) is not equal to one, in fact Tr(A4) = d, the proof of 



which is in Appendix A. 6. 3 



If we repeat the example from Section 6J2 using preparation by projective measure- 
ments with a CNOT gate and a maximally entangled input. We calculate the bilinear 
process matrix to be 

/ En E 12 
E21 E22 
[0] [0] 

V [o] [o] 



M 



[0] 
[0] 

-^22 
E\2 



[0] \ 

[0] 
E21 
En ) 



where Eij = and [0] is the matrix with all entries zero. M. has two non-zero 

eigenvalues Ai = A2 = 1, and Tr(A^) = 2. 

Since Ai is a large matrix, d 3 x d 3 performing tomography is a far more excruciating 
task, as we now require y(<i 2 + 1) input states to completely characterize M. instead of 
the d 2 required for linear process tomography. We did not get time to fully investigate the 
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actually procedure of performing bilinear tomography, and investigating this is a possible 
avenue for future research. It would be interesting to see how initial correlations can 
effect a bilinear assignment map, and under what conditions the evolution would always 
be completely positive. We suspect that similarly a condition for complete positivity 
would be requiring that the bilinear process matrix M. is positive. 
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In the last chapter we provided examples of idealized SQPT experiments where state 
preparation resulted in a non-completely positive process maps in the presence of initial 
correlations. The problem is that in process tomography non-CP results are frequently 
observed due to statistical noise. For example, in the case of photonic qubits this is due to 
the Poissonian count statistics from the spontaneous sources used for creation of photons. 

Often when an experimenter identifies a non-CP process map in a tomography exper- 
iment it is assumed to be due to noise, and optimization techniques such as maximum 
likelihood estimates are employed to map it to the closest physical process map, i.e. a 
completely-positive one. However, if we accept that non-CP maps can occur for reasons 
other than statistical noise, we need some way to distinguish between non-physical re- 
sults due to the noise, and those which arise legitimately due to initial correlations. In 
this chapter we investigate this issue for an optical implementation of quantum comput- 
ing where the statistical noise is Poissonian. In the following chapter we will outline an 
optical experiment to demonstrate the results of this thesis. 
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7.1 Statistical Noise in State Tomography 

We start by considering state reconstruction of a single qubit. This will allow us to build 
an intuitive picture for the statistical noise present in state reconstruction before we move 
to the two-qubit case needed for describing a single qubit process matrix. 



As mentioned in Section |5.3.1[ due to the spontaneous nature of photon sources the 
distribution in recorded counts in an experiment will be Poissonian. If we are measuring 
a single qubit with density matrix p, the expected number of counts for Af measurements 
of an operator M m is 

n m = J\fTr[pM m ]. (7.1) 

The density matrix p is reconstructed from the expected counts using the dual basis, 
{D m }, of the measurement operators, and expected probability coefficients p m = n m /Af 
recovered from the relative frequencies n m . The reconstruction is given by, 



P 



J2Pm D ^ ( 7 - 2 ) 



The influence of statistical noise means that for any given tomography experiment, 
the measured count data will be a set of random variables obeying Poissonian distri- 
butions P(A) with the distribution parameter given by the expected number of counts 



from Eqn (7.1). For example, the distribution of the measured counts for measurement 
operator M m is n e m ~ P(n m ). 

We will assume that our experiments run for sufficiently long to approximate this with 
a Gaussian distribution N(p, a 2 ) with mean and variance (p m , cr^) = (n m , n m ) [35]. Hence 
the distributions of the observed counts are given by 

n e m ~N(n m ,n m ). (7.3) 

This approximation is very accurate for large values of n m , and hence is most valid 
when dealing with large M . However, its accuracy will be reduced for states p which are 
orthogonal or close to orthogonal to one of our measurement operators, Tr (pM m ) ~ 0), 
as in these cases the expected count number will be very low. By making this approxi- 
mation we can use some useful properties of the Gaussian distribution. For a normally 
distributed variable X ~ N(p,a 2 ) and scalar a, aX is normally distributed with mean 
p! = ap and variance a 2 ' = a 2 a 2 . Also, for a set of normally distributed random variables 
Xi ~ N(pi,a 2 ), by the Central Limit Theorem [22] the sum of the variables is normally 
distributed with mean p = YliVi and variance a 2 = J2i °f • 
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Hence we have that the measured values of the probability coefficients from Eqn (7.2), 
Pm = n mlN ■, are normally distributed variables 

Pi ~ Nfe^) (7.4) 



N (7.5) 



N ' M 2 

Pm 
Jf 

where p m = n m /Af is the expected value of p e m . 

Combining these we obtain our distribution for density matrices reconstructed via 
state tomography. The measured density matrices p e will be distributed according to 



J2^(p m ,jf)D m (7.6) 



which is a multidimensional sum of Gaussian distributions, as the dual matrices D m define 
different directions on the Bloch sphere. 

If we examine the distribution in a specific dimension, for example X = a±, then 



Tr[Xp e ] ~ ^Tr[XD m ]N(p m , 



Pm 

M 



N (p m Tr[XD m \, ^Tr[XD r , 

m 

N \J2 Pm Tr[XD m ], J2 ^Tr[X,D, 

\ m m 

N (tt[X P ], -^^^TrfX^] 2 ] , (7.7) 



which is a Gaussian distribution with mean /x = Tr(Xp) and variance a 2 = ^ m p m Tr[XD m ] 2 /jV. 
So the variance of the distribution scales as 1/jV. This holds for all dimensions. 



In Fig. (7.1) we plot the reconstructed density matrix for a p = \D)(D\, for 1000 
repeats of the reconstruction process. This is done for M = 10 2 , 10 3 , and 10 4 . Here we 
can see that the distribution is indeed centered on p, and has the appearance of a three 
dimensional Gaussian ball. We can see that the variance decreases with M and that 
roughly half the data represents unphysical states, as they lie beyond the surface of the 
Bloch sphere. 



7.2 Statistical Noise in Process Tomography 

Now that we have an understanding of the single qubit case, we move to a process map for 
a single qubit. By the Jamiolkowski isomorphism reconstructing a process matrix for a 
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|H> |H> |H) 




(a)AA=10 2 (b)AA=10 3 (c)AA=10 4 

Figure 7.1: Distribution of reconstructed states due to Poissonian noise for p = \D)(D 



single qubit is equivalent to reconstructing a two-qubit density matrix. A two-qubit state 
is characterized by 12 independent parameters, and the reconstruction procedure requires 
16 independent measurements, instead of the four needed for a single qubit. We can follow 
the same procedure as the single qubit case and reconstruct a probability distribution of 
output states, though we no longer have a convenient visualization. 



For the reconstruction of a process matrix we can generalize Eqn. (7.1). The 
expected number of counts when measuring input state pj with measurement operator 
M m is 



A/"Tr A [{pj ® M m )A £ ] (7.8) 
The excepted reconstruction for the process matrix is then 

he = ^p jm D* ® D m (7.9) 



where pj m = rij m /Af is the expected probability coefficient for D* eg) D m , and {Dj} and 
{D m } are the dual bases to the input basis {pj} and measurement basis {M m } respec- 
tively. We note that if our measurement and input bases are the same, so too are the 
corresponding dual bases. 

It is then trivial to extend the distribution of reconstructed density matrices to process 
matrices. The distribution of measured process matrices A| is 

A| ~ J2 N {Pi™ ) D *j ® D m- (7-10) 

jm 

which is a 16-dimensional sum of Gaussian distributions. 
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If we examine the distribution in a particular dimension, say X ®Y for example, then 



analogous to Eqn. (7.7) it is given by the Gaussian distribution 

Tr[(X®Y)A|] ~N(^,(T 2 ) , 

with mean p = Tr[(X ®Y)A g ] and variance a 2 = ± J2 jm TrtX/^fTrfFAn] 2 - This result 
holds for all dimensions, hence like the case of a single-qubit state, the variances in the 
distribution of a reconstructed the process matrix scales as 1/M in all dimensions. 



7.2.1 Non-CP Process Matrices 

The reconstruction distribution for a non-CP process matrix will be identical to the CP 
case, with respect to the appropriate mean and variance. The only difference now is 
that the distribution will be centered outside the physical subspace. Since the variance 
of the distribution scales as this gives us a possible way of distinguishing between 

non-CP cases due to noise, and those arising legitimately. If Af is sufficiently large, then 
an overwhelming majority of the distribution from a non-CP process matrix will also be 
non-CP. However in the case due to noise, a larger portion of the distribution will still 
correspond to CP processes. With this assumption, we should be able to employ standard 
statistical hypothesis testing to determine, within a confidence bound proportional to A/", 
the probability of a measured non-CP result being due to noise. 

We can illustrate this idea using a single qubit density matrix as an example. Consider 
two single-qubit density matrices pA and ps, with Bloch vectors cxa = (l/v^2, 0, l/\/2) and 
Ob = (1-2/ V%, 0, 1.2/v^2) respectively. These two vectors share a common direction on 
the Bloch sphere, however while pa is a pure state on the surface of the Bloch sphere, pB is 
not a valid density matrix. It has a greater than unit Bloch vector and thus lies outside the 
physical space of the Bloch sphere. We see here that pB is a non-positive density matrix, 
it has a negative eigenvalue of —1/4. We plot the distributions of reconstructed density 



matrices based on these expected states in Fig (7.2). We can see in Fig (7.2(a)) that for 



M = 100 there is a large degree of overlap between the two distributions. However by time 



time M reaches 1000, as in Fig. (7.2(c) ) we can distinguish between the two distributions 
with a reasonable level of accuracy. The case of the process matrix is essentially the same 
however we don't have such a convenient visualization. 
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Figure 7.2: Distinguishing between reconstruction distributions of two density matrices with 
Bloch vectors in the same direction, one positive (black distribution) and one non-positive (blue 
distribution) . 



7.3 Summary and Future Work 

We require a method of testing a measured non-CP process matrix to assign probabilities 
for the deviation from complete positivity being due to statistical noise and being due 
to initial correlations. As we previously mentioned, the distribution of process matri- 
ces reconstructed by process tomography (both CP and non-CP) is a multi-dimensional 
Gaussian distribution. In each dimension however the variance is inversely proportional 
to A/", so by greatly increasing M we can increase the distinction between the two causes 
for a non-CP result. However, in certain situations, such as a non-CP process map lying 
very close to the CP boundary, increasing M to the required level may not be realistically 
achievable. 

A more rigorous formulation of a statistical model to allow for non-CP process ma- 
trices in tomography optimization would be an entire study in itself. This would be an 
important direction for future work in the area of dealing with non-CP results in process 
tomography. A possible direction this research could take is to consider more advanced 
sampling schemes for our reconstruction data, such as sequential analysis or adaptive 
schemes. 

In sequential analysis the sample size M is not fixed in advance. Instead data is 
evaluated as it is accumulated. We start with a certain sample size A/i, and perform 
tomography measurements, then following reconstruction we repeat with process with 
another sample A/2- This count data is added to N\ and reconstruction is repeated with 
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the cumulative data. This process repeats until some predefined termination rule when we 
can reconstruct a significant result us satisfied. A possible termination rule is reaching a 
high confidence level in identifying the measured process matrix as a true non-CP result, 
or as an artifact of statistical noise. 

In an adaptive scheme we would change our measurement procedure over the course 
of the experiment. For example, if we combined an adaptive measurement scheme with 
sequential analysis, each time we collect a new set of count data we could change our 
measurement basis. This would make the reconstruction algorithm more complicated, 
however it has the potential to increase the accuracy of our results. 
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8 

Experimental Implementation With Linear 

Optics 



The major part of this thesis has been focused on the mathematical description of open 
quantum systems, and their evolution in the presence of initial correlations. In this 
chapter we move into a physical picture where we impart to the reader a sense of how 
such abstract concepts might be demonstrated in an actual experiment. We will focus on 
an implementation using linear optical quantum computing (LOQC) with photonic qubits. 

LOQC is one of several major architectures currently used to implement quantum 
information techniques [34J. It is sufficiently advanced to perform quantum gates, it has 
well characterized sources of noise, and well developed measurement and tomography 
techniques. However, a thorough description of LOQC is far beyond the scope of this 
thesis. Instead we will give a very brief introduction to some of the key concepts involved 
in translating abstract mathematical entities into actual physical devices. Our goal is to 
give an outline for an experiment designed to implement some ideas proposed in this thesis. 



Specifically, we are concerned with implementing the example shown in Section 6.2 which 
highlighted how different state preparation techniques could lead to to CP or non-CP 
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results when used in SQPT. 



8.1 Photonic Qubits 

In LOQC there are two main methods of encoding a photon as a qubit. The first, which 
we mentioned in Section 2.3\ is to use the polarization state of the photon. Here the logical 



states |0) and |1) are taken to correspond to the horizontal and vertical polarization states 
respectively. The second method, known as the dual rail representation, is to represent a 
qubit by a photon in one of two separate optical modes. The logical states |0) and |1) are 
taken to correspond the photon being in the first mode |1,0), or the second mode |0, 1) 
respectively. We should note that polarization encoding can be considered as a specific 
type of dual rail encoding. The implementation for performing one and two qubit gates 
that we will examine is based principally on polarization encoding. 

To generate photonic qubits, single photon sources are required. Currently this is 
done by a process known as spontaneous parametric down conversion (SPDC) [lj. This 
process involves pumping a nonlinear optical medium with light of frequency ujq. With 
some small probability a pump photon will spontaneously down convert into two daughter 
photons, each of frequency ujq/2. The careful implementation of filtering allows us to use 
these daughter photons in an experiment. However, since the down conversion process 
is spontaneous and photon detection is destructive, to use both the daughter photons 
an experimenter will usually employ post-selection techniques. This means that only 
instances when two photons are detected at the end at the end of the experiment will be 
considered valid. 



8.2 Performing Single Qubit Unitary Operations 

LOQC has everything we need to implement arbitrary single qubit gates. This is done 
through the use of wave plates, phase delays, beam splitters, and polarizing beam splitters. 

For polarization-encoded photons, single qubit operations can be implemented using 
combinations of wave plates. These are optical devices made from birefringent crystal 
which induce a phase shift between orthogonal polarization components of a light beam. 
In practice we only need half- wave plates (HWP) and quarter- wave plates (QWP), which 
induce a phase difference of <p = it and <p = it/ 2 respectively [38J. The unitary operators 
corresponding to HWP and QWP with the axis of the induced phases shift rotated an 
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angle 9 from horizontal are given by 



Uhwp(9h) 



cos 28h sin 29^ 
sin 26h — cos 29^ 



U, 



qwp\ u q 



V2 



l + i cos 26 q 
i sin 29 n 



i sin 29 q 
1 — i cos 29 n 



From here any single qubit unitary gate can be implemented on a polarization encoded 
photon using a combination of appropriately aligned HWPs and QWPs. 



8.3 Two-Qubit Entangling Gates 



To actually perform quantum computation we required a two-qubit gate which can induce 
(or remove) entanglement between a pair of photons. Such a gate is called an entangling 
gate. Implementing a two-qubit entangling gate is more challenging than a single-qubit 
gate, as inducing interactions between photons is quite difficult. This issue has been 
circumvented by using a process called measurement induced non-linearity [39]. In effect, 
this simulates a nonlinear interaction with a measurement. The trade off here is that the 
process is non-deterministic, which means there is only a certain probability of the gate 
being successful, and success cannot be determined without some form of measurement. 

Conceptually, the simplest such gate is the controlled-Z (CZ) gate. This gate induces 
a phase shift of tc on the |11) term of a two-qubit state. Its action is given by the unitary 
matrix 

" 1 
10 
10 
-1 

Recently three groups [1DJ HH H2] independently arrived at similar designs for a particu- 
larly simple implementation of a CZ gate in LOQC. The implementation by Langford et 
al. (ID] requires only three partially polarizing beam splitters and two half-wave plates. 



U, 



cz 



The design of their gate is shown in Fig. 8.1 



CZ gates in combination with single-qubit gates can reproduce any two-qubit gate pQ . 
To implement a CNOT gate we require the application of two Hadamard gates, one 
applied to the target qubit before the CZ gate, and one after. The unitary matrix for the 
Hadamard gate is given by 



H 



1 

71 



i i 
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7 * 

1/2-wave plate [] swap A/2 @ 45° [XI P artial PBS (>lh=1/3, r|v=1) 
1/4 -wave plate \Z\ PBS -^{HZM* polarization analyzer 

/ beam-splitter (reflectivity n=1/3) single-mode fibre launcher 

Figure 8.1: Linear optical implementation of a non-deterministic Controlled-Z gate gate by 
Langford et al. [30]. Figure reproduced with permission. 

which up to a global phase shift of tt is simply a HWP with 9 set to 22.5° from horizontal. 
With the necessary tools for producing qubits, performing a CNOT gate, and determining 
the output states with SQPT, all that is left is the initial state preparation. 



8.4 State Preparation in LOQC 

In comparison to the theory presented in Chapter [5j the only novel factor introduced 
in an experimental implementation of SQPT is the method of state preparation for the 



tomography input states. In Section 6.1.2 we proposed two different methods of prepa- 
ration which involved projective measurements, one using only measurements, and the 
other using a single measurement followed by unitary rotations. We will now recast these 
procedures in terms of optical components. We are assuming that our experiment consists 
of two photonic qubits which can be prepared in a correlated initial state 70. We will call 
the first qubit the system, denoted S, and the second qubit the environment, denoted E. 
Our preparation procedures are only applied to the system qubit S while the environment 
qubit is left alone. 
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8.4.1 Method I: State Preparation by Measurement and Rota- 
tions 



First we will consider the preparation map described in Section 6.1.2 Recall that it con- 
sisted of performing a single projective measurement to a fixed pure state, and then using 
unitary rotations to transform the post-measurement state to the required inputs states 
for performing SQPT. We wish to prepare qubit S into the input states \H) , \ V) , \D), 
and \R). 

Our first problem with implementing this procedure in LOQC is that we cannot truly 
perform projective measurements. This is because measurements in optics rely on photon 
detection, which destroys the photon being measured. One way to simulate the effect of a 
projective measurement is to use a polarization filter and to post-select photons that get 
transmitted as verified by a final detection. This filter can be rotated to only transmit 
light of any chosen plane-polarized state, the unitary matrix for a polarization filter with 
its polarization axis set at an 6 from horizontal is given by 



1 + cos 26 p sin 26 p 
sin 26 p 1 — cos 26 p 



If we align a polarization filter to transmit only \H) (6 = 0), then we can use a HWP 
and QWP to give our desired input states (relative to a total phase shift experienced by 



the qubit). The preparation apparatus for this procedure is depicted in Fig. (8.2). For 
Method I the polarizer is fixed at P(0). The alignment of the HWP and QWP required to 
give our required input states, and the induced phase shift of the wave plates, are shown 



in Table 8.4. 1| We note that when using this preparation procedure for SQPT we always 



expect to reconstruct a CP process matrix. 



8.4.2 Method II: State Preparation by Measurement 



Now we consider the second preparation procedure, which we outlined in Section 6.1.3 
This procedure described preparing our four input states {\H) ,\V) ,\D) ,\R)} by only 
using projective measurements. As with Method I, since we cannot perform a true pro- 
jective measurement, we simulate the effect of one using a linear polarizer. However, we 
now have an additional problem as we cannot directly produce elliptically polarized light 
in such a manner. We propose then that we overlook this problem and prepare the state 
\R) in the same manner as in Method I. For the vertical and diagonal states however, we 
can change the angle of the linear polarizer instead of using the HWP and QWP as in 
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Source P HWP QWP 




Figure 8.2: Apparatus for performing state preparation in linear optics for Methods I and 
II. Settings for the angles of optical elements to prepare the input states H,V,D,R are shown in 



Table (8.4.1). 





Method I 






Method II 








9 P 


o h 


o q 


z-shift 


P 


o h 


o q 


z-shift 


\H) 


0° 


0° 


0° 


135° 


0° 


0° 


0° 


0° 


\V) 


0° 


45° 


0° 


45° 


90° 


0° 


0° 


225° 


\D) 


0° 


22.5° 


45° 


135° 


45° 


22.5° 


22.5° 


135° 


\R) 


0° 


45° 


-45° 


0° 


0° 


45° 


-45° 


0° 



Table 8.1: Settings for optical implements to prepare required input states for SQPT for 
methods I and II. Angles are rotations of optical axes from horizontal for a linear polarizer (0 P ), 
HWP 8h, and QWP 8 q . z-shift is the relative phase shift experienced by qubit S compared to 
qubit E. 



method 1. As we will show below, this revised preparation procedure can still result in a 
non-CP process matrix. 



We will use the same preparation apparatus as Method I as depicted in Fig. (8.2). 
However, now the angles on the optical implements will be different. The settings for the 



polarizer, HWP, QWP and the total induced phase shift are also shown in Table 8.4.1 We 
make note that in the case of producing diagonally polarized light, the angle settings for 
the HWP and QWP were not zero, this is because for such an alignment, the combination 
of the wave plates would rotated the diagonally polarized light transmitted by the polarizer 
to anti-diagonally polarized light. We could just remove the HWP and QWP from the 
setup. However they are still needed to produce our right-circularly polarized state. In 
addition using this setup gives us a constant state preparation apparatus consisting of a 
polarizer, HWP, and QWP. 

As previously mentioned we cannot perform a projective measurement for the state 



8.5 Outline of an Optical Demonstration of Non-CP Maps 



73 



\R), and so we produce \R) by rotating \H). We will now revise our SQPT example from 
Section [6^21 to accommodate this. 

We suppose our system and environment, both photonic qubits, are initially in the 
maximally entangled state \<f>) = ( 1 00) + |ll))/v2. The interaction between the two 
photons is given by a CNOT gate with the system photon the target qubit. If we perform 
SQPT with this hybrid preparation procedure the resultant process map is given by 





2 





-1 -i 


2 


1 











l + i 


2 


-1 + z 





2 







2 


1-i 









which has eigenvalues {2.039, —1.039, 0.863, 0.137}. Hence we still obtain a non-CP result 
with this method. 



8.5 Outline of an Optical Demonstration of Non-CP 
Maps 

The previous sections illustrate all the components needed to perform an experiment to 
demonstrate non-CP maps, which could be implemented with currently available tech- 



niques in LOQC. A simplistic schematic for the experiment is shown in Fig. (8.3). The 
details on preparing the initially correlated state, and performing state tomography at 
the output are to be determined according to the planned experiment set-up. 

Following the preparation of an initially correlated joint system state 70, the exper- 
imenter will perform SQPT of the CNOT gate, which is constructed using the CZ and 
Hadamard gates. The input states for SQPT should be prepared using both Methods I 
and II as perviously outlined, and the reconstructed process matrices should be compared. 
For the initially correlated states we suggest the experimenter first prepare the maximally 
entangled Bell state 70 = \<j))(4>\. Following this the experimenter should prepare separable 
states, for example 70 = \ {pu ® Pa + Pd ® Pv), and 70 = \ {pu ® Ph + Pv ® Pd)- 

8.5.1 Analysis of Results 

Using the mentioned initial states 70, our expected results are a CP process matrix for 
Method I, and a non-CP process matrix for Method II. Problems may arise in SQPT 
because we can no longer use maximum likelihood tomography since we are explicitly 
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Figure 8.3: Design of an experiment to demonstrate non-CP maps with currently available 
technology in linear optics. 

looking for the possibility of a non-CP process map. To reduce the influence of statistical 
noise in the count statistics we suggest running the experiment for much longer than 
usually done so as to minimize the effect of the Poissonian noise on the count data. Based 
on the statistical methods introduced in Chapter [7] we should increase A/", the number of 
copies of a state used in the tomographic reconstruction, to as large as feasibly possible. 
This will increase the ability to determine the probability that a measured process matrix 
arose from our expected process matrix due to statistical noise. 
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Conclusion 



9.1 Summary 

In this thesis we examined the evolution of open quantum systems in the presence of initial 
correlations in relation to non-completely positive maps. The study of open quantum 
systems is a crucial subject in quantum information science, as all real world quantum 
devices will experience some degree of interaction with their surroundings. In the case 
where we are not interested in the explicit time evolution of such systems we describe 
their evolution by quantum operations. These are CPTP maps on density matrices. In 
the early chapters we examined the theory for CPTP maps. 

In Chapter [4] we introduced non-CP maps. These can occur in quantum systems when 
the state of the system and its environment are initially correlated. In this chapter we 
investigated a new scheme of classifying the correlations of quantum systems by using 
a quantity known as quantum discord. We found that the proposed new definition of 
a classically correlated state being one with zero quantum discord is inconsistent. We 
showed, by a counter example, that not only is the quantity of quantum discord not 
symmetric, but if the state of a joint system has zero quantum discord, then the discord 
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of the state with the roles of the two systems interchanged in generally need not be zero. 

If we define correlations as quantum when a state has quantum discord greater than 
zero, this can lead to a situation where a system is both quantum and classically correlated 
- two situations which intuitively should be mutually exclusive. This result implies that 
when using the proposed definition of classifying correlations based on quantum discord, 
correlations in quantum mechanics are directional. We also investigated a result which 
stated that the evolution of initial states with only classical correlations would always lead 
be CP. In light of the found asymmetry we a led to conjecture that this theorem would 
not hold if we reversed the roles of the systems in our definition of quantum discord. 

In Chapter [5] and Chapter [6] our research moved into the direction of quantum pro- 
cess tomography. We introduced the two main schemes of implementing this important 
technique for completely characterizing an unknown physical process: SQPT and AAPT. 
We investigated the effect of state preparation on these procedures, examining methods 
proposed in the literature, and also proposed our own, which has number of advantages. 
We investigated how these procedures influenced the nature of the reconstructed process 
matrix, either CP or non-CP, an question which has not been asked before. 

We found that when using AAPT to characterize an unknown quantum operation, the 
reconstructed map would always be CP regardless of initial correlations. This is because 
only a single input state was required, effectively truncating the influence of the environ- 
ment on the system by fixing it in a single state in accordance with the standard theory of 
CP maps. In the case of SQPT however, the nature of the preparation method was found 
to have a large influence on the nature of the reconstructed map. We investigated to spe- 
cific examples which could be implemented in linear optics, a preparation procedure from 
the literature which used only projective measurement, and a procedure of our own design 
which used a single projective measurement followed by unitary rotations. We found that 
our proposed method would always lead to the reconstruction of a CP map, however the 
procedure using only measurements was found to lead to non-CP evolution. It is known 
that the preparation by measurement is a bilinear quantum operation, as apposed to 
linear operation. We proposed a novel method of implementing the bilinear preparation 
procedure as a linear process, which does not affect the issue of non-CP dynamics. 

At the end of Chapter [6] we investigated a proposed procedure for describing bilinear 
process matrices. We noticed that one of the properties of the bilinear process matrix 
asserted in the literature was incorrect. Namely, the author stated that a bilinear process 
matrix has a trace equal to one, when in fact the trace is equal to the dimension of the 
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system undergoing evolution. 

An important issue which arose in our investigation of process tomography was the 
effect of statistical noise on the reconstructed process matrices. We examined the case 
of quantum optics where the spontaneous nature of photon emission leads to Poissonian 
distributed count statistics. This statistical noise could lead to the reconstruction of non- 
CP process matrices when we expect a CP result. This presents a challenge of having to 
distinguish between non-CP results due to noise, and true non-CP reconstructions in any 
actual tomography experiment. 

In Chapter [7] we formulated a statistical distribution for reconstructed process ma- 
trices. This was given by a multidimensional sum of Gaussian distributions, each with 
variance inversely proportional to the number of copies of our states used to perform 
tomography. We proposed that by increasing the number of copies one could distinguish 
between true non-CP results and ones due to noise, as the distribution for true non-CP 
results would largely lie outside CP process space, when in the noise case the majority of 
the distribution would correspond to CP processes. 

Finally, in Chapter [8] we brought together many of the ideas proposed in the earlier 
chapters to design an experiment which could demonstrate non-CP maps. This was done 
using currently available equipment from linear optics, and demonstrated how different 
state preparation procedures in SQPT could greatly influence the identified quantum pro- 
cess. Our two proposed preparation procedures both consisted of three optical elements, 
a linear polarizer followed by a HWP and QWP. The only difference between the schemes 
was the orientation of the optic axes of each device. The successful implementation of 
this experiment would illustrate how non-CP makes can arise in real world situations, and 
also show how a seemingly minor change in preparation can have a large impact on our 
final tomographic output 

9.2 Directions for Future Investigation 

Our research left several avenues for further study into non-CP quantum processes. An 
obvious step would be to actually implement our proposed experiment from Chapter [8j 
In addition a more rigorous study of the statistical distribution of reconstructed density 
matrices would be beneficial in helping use distinguish between true non-CP results, and 
ones due to statistical noise. Some possible directions for this investigation include de- 
veloping a tomography protocol which incorporates sequential analysis and or adaptive 
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schemes. One possible example for this would be an implementation of tomography where 
the number of copies prepared is not fixed, but we can resample as required to increase 
accuracy, changing the direction of our measurements to minimize the variance of the 
resulting statistical distribution. Finally, our study of quantum discord in Chapter [4] 
illustrated inconsistencies with the proposed method for classifying correlations of quan- 
tum systems. While the approach appears to have merit, the lack of symmetry poses an 
interesting question about the symmetry of system correlations. 




Mathematical Proofs 



A.l Proofs for Chapter 3 



A. 1.1 Process Matrix Equations 



We wish to prove that for defined by A^ = . (g> S(Eij), the quantum operation 
£(p) is given by Equation (3.5). Note that we can expand any density matrix in terms of 
Eij as p = PijEij. Hence p T = PijEji- Then the right hand side of Equation (3.5) 
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Tr A (A(p T ®I)) = Tr A 



is, 

E ij ® £{ E ij) I I 5^ PmnEnm ® I 
ij / \mn / 

= PmnT^A [{E tj E nm ) <g> 

ij rrtn 

= ^ ^ Pmn£(Ejj) Tr [£'^-£' nm ] 

ij mn 

= EZ)^W&[l<>0'|n>Hl 

ij mn 

^ ^ ^ ^ £(.PmnEij^)5j n 6j m 
ij mn 

ij 

= £{p) 

A. 1.2 Process Matrix to Kraus Representation 

We wish to construct a Kraus representation for a quantum operation £ described by a 
process matrix A^. To do this we first prove that for a basis E^ = we can express 

evolution by the quantum operation £ as, 

£(p) = J2(^£)rr>,ss>E? r , P E ss ,. (A.l) 
rr' ss' 

This is called a canonical Kraus representation. 

Our first step is to expand p in terms of E^, as p = PijEij- Hence 

Ej r ,pE ss i = pijE r i r EijE ss / = p rs E r i s i 

ij 

Next, from the definition of A we have 

{he)rr',ss> = ^2[Eij®£(Eij)] rr ',8s' 

ij 

^ ^ (Ejj^ rs £(Eij^ r ' s i 

ij 

^ ^ £ (.Eij ) r > s i 5i r 5j s 

ij 

£{^E rs ^ r i s t 
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In addition, we notice that we can expand out £(p) as 

£(P) = ^ £ (PrsErs) =^ £{p)r's> = ^ £{p rs E ra ) r 
rs rs 

Finally we put these three facts together to get 

y^ J (^£)rr',ss'Ej r ,pE ss i = Prs£ (E rs ) r > s > E r i , 

rr' ss' rr' ss' 

= £{p rs E rs ) r i s iE r i, 

rr' ss' 

= S(p), 



so equation (A.l) does indeed describe evolution by £. 



Now we construct the Kraus operators K n . First we collect the double indices in 



Equation (A.l) to single indices, so £(p) = ^2 mn {As)mnE m pEl l . Next, since the process 



tK e i)n-En 



matrix is positive, we can find its spectral decomposition As = £\ Aj |ej)(ej|. Hence, 

J2(^s)mnE^pE n = J2J2 X ^^nE m pE n 

mn i mn 

i \ m / \ n 

i 

where Ki = Ylm \^i( e i)mE^, and (e$) m is the m th element of the i th eigenvector of Ag. 
Finally, we notice that J2 m V / ^i( e i)m-E'm can ^ e rewritten as y/Ximat |e») . 

Hence we have shown that for £ defined by a process map As, we can construct a 
Kraus representation by forming the Kraus operators Ki = y/XiTnat |ej), where Aj and 
\ei) are the eigenvalues, and eigenvectors of A^ respectively. 

A. 1.3 Superoperator Definition 



We will show that our definition of the superoperator in Eqn (3.7) is consistent with the 



evolution described by Eqn. (3.6). First we note that the set {Eij} is an orthonormal 



basis with respect to the Hilbert-Schmidt inner product. To check 

(EylEu) = Tr(4.£ fe 

= Tr(\j){i\k)(l\) 
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Hence, any density matrix can be written p = J2ijPijEij, or i n vectorized form \p) 
YliiPij So the evolution of a vectorized matrix \p) is given by 



*e\p) = ^{^PiAE^ 
= ^Pij^slEij) 



13 

ij kl 

= ^2pij^2\ S ( E kl))SkiSlj 

ij kl 

= », sii-:,,)} 

= Ysisip^)) 

ij 

= \£{p)). 

A. 1.4 Product states with fixed environment always give com- 
pletely positive dynamics 

Assume our joint system is in a simply separable state 7 = p <8> t, where p is an arbitrary 
state of system A and r is a fixed state of system B, independent of the state of system 

A. We will assume system AB evolves via an arbitrary unitary operation U . To prove 
complete positivity we will find a Kraus representation for the quantum operation £ given 
by =Tr s [U{p®r)W]. 

Since r > we can take its square root, so the initial state is 

7 = P ® VtVt = (I <S> y/r)(p ® I) (I ® v 7 ^)- 
We can then expand out the middle I as £V where {\i)} is an orthonormal basis for 

B. Hence 7 = ]T\(I ® v^)(P ® |«> <«j) (K ® v 7 ^)- 
Our evolution is then given by 

£{p) = Tr s [[/ 7 C/t] (A.2) 

= J^Ttb [U(I®yfr)(p®\i)(i\)(I®^)tf] (A.3) 

i 

= Y(j\ B U(^\i)) B p((i\V^) B rf\j) B (A.4) 
= ^ K tJ pKj r (A.5) 
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where we are using the vector description of the partial trace from section (2.8), and the 
superscript B denotes an operator acting on system B. 

So we have a Kraus representation given by operators = (j\ B U{y/r \i)) B ■ Now we 
just need to check these operators satisfy the completeness relation. 



E4^. 



Hence £ is CP 



£«*lv 

ij 


^) B tf \j)(j\ B u(V^\i)) B 


(A.6) 


i 


^) B u^u{^\i)) B 


(A.7) 


Em* 

i 


t \i) 


(A.8) 


IaTt(t) 




(A.9) 


I A . 




(A.IO) 
□ . 



A. 2 Proofs for Chapter [4] 

A. 2.1 T>b-.a{p) = implies completely positive evolution 

Consider a system A with its environment given by the system B. We wish to prove 
that if the initial state, pab, of the joint system AB is classically correlated as defined by 
T^B:a{p) = 0, then evolution of the system A will always be CP. 



Let pab satisfy T>b-a{pab) = 0. According to equation (4.12) this means there exists 
a set of orthogonal projectors {IP/} such that 

PAB = J2 U tP^t 



where Tj = Tr A^jPAB^f) and Pj = Tr(ILy pab)- The initial state of system A is 

7] = Tr B (p A B) = E^ n i- 

j 

Hence evolution of rj is given by a unitary operation U on the joint system as follows: 



r)f = Tr B (uiY^PjUj^r^U^j (A.ll) 
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where the partial trace is expressed in terms of an orthonormal basis for state 

vector space of B, and \n) B = I® \n). Since Tj is a positive matrix we can take its square 
root y/fj, and write it as Tj = y/fjl ^frj = ^Jr~j \m)(m\ Jr]. Hence we have 



rjl = ^^^^pj (n\ B U(Ilj y/r]\m)(m\ y/fj)U" f \n) B 

j n,m 

= E w B u ^ \ m ^ B ^ (H ^ B{J] ^ B 

j n,m 

= Y,Y,P3 D nrn U j D n™, (A.13) 



n,m 



where D 3 nm = (n\ U{^f~j \m)) is an operator acting on system A only. 

Now we wish to remove the j dependence from D J nm . We do this noting that D 3 m 



Y^i D nm Sij, and so we can rewrite equation (A.13) as 



V = E E Pi ( E D nJ^ 3 H A U M D nm , (A- 14) 

j n,m \ I / V I / 

where n^H/H,- = H, as IT,- are projectors. 

Now since {Ilj} are orthogonal, SijUj = n^H, Hence we have 

= EE^fE^n/Jn^nJ^n^) (a.is) 

j n,m \ I / \ I / 

j n,m 

= J2K aV Kl (A.17) 

a 

where we collapse (n, m) to a single index a, and the Kraus operators are defined 



K a = K nm = E^Ln* = E (nfUi^HfllL 



Now all that remains is to show the Kraus operators satisfy the completeness relation. 
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a n,m l,k 

= Y<J2 u i(( m \V^ B tf\n)(n\ B u(^ k \m)) B n k 



n,m l.k 



2^n,(<m|^ fl ^tf(0*N) B n fc as ^|n>(n|=I fl 

m l,k n 

n^IIfc (m| ^/riy/rk \m) as U is unitary 
^ n fc (m\ r k \m) as UJl k = Si k U k 



in k 



= ^n fc Tr(r fc ) 

k 

= I A as Tr(r fc ) = 1, and^n fc = I A . 

k 

So the operators {K a } are indeed Kraus operators, hence the evolution is completely 
positive. □ 

A.3 Proofs for Chapter [5 

A. 4 Calculation of the Dual Basis for (H,V,D,R) 

We wish to calculate the dual basis {Z^jfor the tomographically complete measurement 



set {Mi}, where i = H,V, D, R. Starting with Equation 5.2, if we vectorize our matrices 
we have 



(M t \D ) = 5a 
M;) (Mi\Dj) = 5 i} \MJ 

m 

\Dj) = M' 1 \Mj) , 



1 — 11 \ — 1\ j 1 
£|M<)<M<|J \D 3 ) 



where M. = ^2 i \Mi)(Mi\, and M. 1 denotes the Moore-Penrose generalized inverse of 
ME 

1 The generalized inverse for a matrix M is a unique matrix M^ 1 satisfying MM~ X M = M, 
M^MM' 1 = M~ x , {MM~ l ) ] = MM' 1 , {M~ 1 My = M~ X M. In the case where M is an in- 
vertible matrix, its generalized inverse is the usual inverse matrix. 
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For our chosen measurement set, 



M 



\M H )(M H \ + \M V )(M V \ + \M D )(M D \ + \M R ){M R \ 



1 

4 



10 







6 1 -i 
1 + i 2 
1-i 



+ 



2 1-i 1 + i 
Hence taking the inverse gives 

AT 1 = 









1 

1 + i 2 

1 + i 

2 1-i 
6 



+ 



1111 

4 4 4 4 
1111 
4 4 4 4 
1111 
4 4 4 4 



1 1 
4 4 



1 

4 J 



+ 



2 


-1 + i 


-1 - i 





1-i 


6 


2i 


-1-i 


1 + i 


-2i 


6 


-1 + i 





-1 + i 


-1-i 


2 



Now finally, we use this in Dj = mat(M. 1 \Mj)) to get 

1 

Di 



Do 



1 


2 


-1+i 


2 


-1-i 





1 





-1+i 


2 


-1-i 


2 



£>3 



.Da 



1 

1 

-i 

1 



A. 5 SQPT Process Matrix Definition 

We wish to show that the definition of the process matrix in SQPT 



i 1 

4 4 

_ 1 i 

4 4 

1 _ j 

4 4 

i I 

4 4 



is equivalent to the definition in Section 3.2.2 We do this by verifying Equation (3.5) 
with this new definition. 

First we express p in terms of our input states, ie p = J2 k PkPk- Then the RHS of 
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Equation (3.5) becomes 



Trg [A^(I <8> p T )~\ = Tr, 



E^;®^)) (e 



Pkpl ® I 



3 



3 fc 



(jPj)Tr PkDj pk and Dj are hermitian 



j k 



^/.Pk^jpj^jk definition of the -Dj 

j k 



n? 



Since this is true for all p we must have = Z}* ® £(pj) 



□ . 



A.6 Proofs for Chapter 6 



A.6.1 A stochastic preparation map does not effect the state of 
the environment 

If we are preparing our system by a map (5 <g> I) (7) = (8> t. Then by combining 

the spectral and Schmidt decompositions, we can express the initial state 7 as 



Hence, 



7 



E Afc l efc ^ 



efc 



^KPkiPk 3 \i£){jk\ ® HfeXif | 



2^)7 = E A ^ fe ^ 5 (l^>0'feD® l**>0'*l 

= |^)(^|®r. 

=^r = Tr A [(~®Z) 7 ] 

= E A *^«l , *>0*l ,& ^l i *>0*l)] 
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If S is trace preserving then 

Tr[S(|^)<^|)] = Tr[|z^<^|]=5, 

Hence, 



k i 



So state of the environment post preparation is independent of the map S, instead it on/?/ 
depends on the initial state 7. 

A. 6. 2 Bilinear process matrix equations 



We wish to show that our definition of the process matrix in Equation (6.13) is equivalent 



to the original definition in Equation (6.12). We do this by expanding out our expression 



into the index notation of the original equation. So, 

A^U' = Tr([f/M 7 oW^L 

= ^2{[ U ]rr'[lo]r"s"[U Ji ] s/s ) et 

c 

= EE([ [/ ]-') ea ([7o] ) -» s ») Q/3 (^ t ]^) /3e 

e a,(3 

= E Ure,r' a{lo)r" f3,s" e{U^) s' p,st, 
a,P,e 

which is the original definition of Ai 



□ . 



A. 6. 3 Trace of bilinear process matrix is d 

Let M be a bi-linear process map given by an initial state 70 and unitary U for a d 2 - 
dimensional joint system AB. Hence M is given by 

M r ( ^, v = l>([^] rr 47o]r» s »[^]L)- 
We now compute the trace of M. 

Tr(m) = £ M^vW 



^(P\rAlv\r»AU\\ r ^ 



Tr 
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But since U is unitary 



Id* = U ] U 



il 



X 



x 



X 



\ 



V 



X 



x EMM** J 

where x are cross terms [^^[t^M- Hence we have 



'J 



Hence 



Tr(m) = Tr 



Oj r"r" 



d ^Tr([7o] r 'V; 



dTr (70) 
d. □ 
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